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OPTICAL  POWER  SPECTRAL  ANALYSIS 
FOR  MACHINE-READABLE  FACTOR  MAPS 


I.  INTRODUCTION 


A  preliminary  investigation  of  the  application  of  Optical  Power  Spectral  Analysis 
(OPSA)  approaches  to  machine-readable  factor  maps  for  use  in  terrain  information 
systems  is  presented.  The  objective  is  to  offer  a  more  compact,  simpler,  and  cheaper  data 
base  storage  and  retrieval  system  than  would  be  available  by  all-digital  technology.  In 
the  basic  concept  of  the  systems,  OPSA  is  used  to  decode  data  from  reduced  factor 
maps  that  have  been  encoded  by  spatial  frequency  and  direction  of  optical  gratings. 


Alternative  concepts  to  high-priced,  large  volume,  and  relatively  slow  digital  mass 
storage  systems  for  operation  in  the  field  are  few.  For  this  reason,  a  preliminary  study 
of  the  application  of  OPSA  was  conducted  in  depth.  The  basic  concept  of  optical  coding 
and  OPSA  decoding  of  terrain  factor  maps  is  presented  in  Chapter  2.  Chapter  3  contains 
an  analytical  characterization  of  the  OPSA  sampling  system  and  its  use  in  detection  of 
Ronchi  gratings.  Chapter  4  presents  a  digital  simulation  program  and  simulation  results 
from  its  use  under  a  range  of  operational  parameters.  The  preliminary  design  of  an  optical 
subsystem  for  detection  of  grating  maps  is  developed  in  Chapter  5  to  show  how  the 
interaction  of  parameters  forces  the  designer  to  make  trade-offs  for  realistic  optical 
system  specification.  In  Chapter  6  the  technology  for  producing  grating  maps  is  reviewed 
and  approaches  for  extending  capabilities  to  produce  gratings  with  increased  storage 
capacity  are  outlined.  Chapter  7  outlines  more  general  system  concepts  which  integrate 
optical,  electronic,  digital,  and  storage  and  retrieval  subsystems  and  an  example  of  a 
grating  detection  system  is  described  to  illustrate  subsystem  interactions.  Chapter  8 
presents  the  conclusions. 


II.  BASIC  CONCEPT 


If  we  insert  a  uniform  grating  in  a  collimated  coherent  light  beam,  as  in  figure  1. 
at  the  hack  focal  plane  of  the  transform  lens,  we  will  detect  a  number  of  collinear  light 
spots.  The  spacing  of  the  light  spots  will  be  determined  by  the  grating  frequency  and  type 
of  grating,  with  the  intensity  of  the  spots  related  to  the  type  of  grating.  The  orientation 
of  the  line  of  spots  will  be  perpendicular  to  the  direction  of  the  linear  grating.  If  we 
change  the  grating,  the  spots  will  change  accordingly  in  the  back  focal  plane  of  the  lens. 
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Note:  This  Uniform  grating  is  illuminated  hy  a  collimated  coherent  light 
beam  that  results  in  collincar  diffraction  spots  in  the  hack  focal  plane 
of  the  transform  lens. 


FIGURfc  1.  A  Uniform  Grating. 
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Knowing  the  type  of  grating  employed  in  the  optical  system,  we  can  determine  without 
error  the  grating  frequency  and  orientation  by  taking  a  few  simple  measurements  on  the 
light  spots. 


Now  consider  a  simplified  attribute  map  as  shown  in  figure  2,  wherein  each 
bounded  area  is  a  two-dimensional  representation  of  a  range  of  one  or  more  properties 
coded  by  an  alphanumeric  character.  Decoding  is  accomplished  by  reading  the  description 
associated  with  the  character  code  in  the  map’s  legend.  Within  each  area  having  the  same 
character  code,  the  properties  are  expected  to  have  the  same  range  of  values,  and  dis¬ 
continuities  occur  only  at  boundary  lines  between  coded  areas.  The  example  of  an 
attribute  map  in  Military  Geographic  Intelligence  (MGI)  is,  of  course,  the  terrain  factor 
map  where  separate  maps  code  surface  configuration,  surface  composition,  vegetation, 
etc. 


*• 

Instead  ot  a  alphanumeric  character  codes,  we  can  use  grating  codes  as  shown  in 
figure  2  with  the  boundaries  retained.  The  code  symbols  are  now  a  two-tuple  (w,  0), 
with  w  representing  an  element  or  the  discrete  set  of  spatial  frequencies,  w  e  w. 
and  0  representing  an  element  of  the  discrete  set  of  grating  directions.  0  e  0. 


The  grating  map  of  figure  2  can  be  a  film  transparency  with  the  grating  lines  having 
a  high  optical  density  and  the  space  between  the  grating  lines  having  a  low  optical  den¬ 
sity.  If  the  optical  grating  map  is  placed  in  the  optical  power  spectrum  sampling  system 
shown  at  the  top  of  figure  3,  that  part  of  the  grating  illuminated  by  the  collimated  co¬ 
herent  light  beam  passing  the  sampling  aperture  will  result  in  spots  of  diffracted  light, 
as  indicated  in  figure  1.  When  a  Recording  Optical  Spectrum  Analyzer  ( ROSA)  detector 
that  has  32  ring  detector  elements  and  32  wedge  detector  elements  with  (see  figure  3) 
each  element  providing  an  electrical  output  signal  is  centered  on  the  optical  axis  in  the 
back  local  plane  of  the  transform  lens,  the  diffracted  light  spots  from  any  grating  will 
result  in  a  set  of  detected  signals  for  those  elements  so  illuminated. 


The  remainder  of  this  report  will  show  that  the  grating  codes  can  be  optimally 
designed  so  that  when  they  are  encoded  in  an  optical  grating  map  and  used  in  a  pro¬ 
perly  designed  optical  system  the  result  is  several  hundred  unique  coding  classes  for 
each  map  sheet,  l.ach  possible  pair  vector  has  a  unique  address  in  a  computer  lookup 
table. 


II 


1 


Note:  The  attribute  map  is  coded  with  alphanumeric  characters  and  the 
grating  map  has  attributes  coded  with  gratings  of  different  spatial 
frequency  and  direction. 

FIGUKt  2.  An  Attribute  Map  and  a  Grating  Map. 
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III.  ANALYTICAL  CHARACTERIZATION 


In  this  chapter  the  analytical  characterization  of  ROSA  and  the  theoretical  aspects 
of  OPSA  as  applied  to  optical  gratings  and  their  detection  will  be  discussed.  This  material 
will  serve  as  a  background  for  subsequent  discussions  involving  a  computer  simulation 
program  that  will  demonstrate  grating  discrimination  potential  with  ROSA  systems  and 
involving  the  design  of  grating  detection  systems.  The  ROSA  system  hardware  has  been 
previously  reported,1  and  the  system’s  use  for  automated  pattern  recognition  research 
is  documented  in  several  publications.2  -1-4 


The  approach  used  here  should  not  be  considered  rigorous,  but  rather  more 
intuitive  so  as  to  provide  a  reader  without  a  background  in  OPSA  with  enough  founda¬ 
tion  to  understand  the  relation  of  the  design  variables  and  the  nature  of  the  trade-offs 
to  achieve  a  useful  operational  system. 


A.  ROSA  Characterization  •  The  three  ROSA  subsystems  are  indicated  in 
figure  4.  Here  we  will  be  concerned  with  the  signals  produced  by  optical  sampling  and 
the  nature  of  the  electronic  signal  processing  in  the  ROSA  system.  The  signal  flow  and 
analytical  characterization  of  the  optica)  and  electronic  signal  processing  within  the 
ROSA  is  summarized  in  figure  5.  We  will  now  develop  the  expressions  shown  in  this 
illustration. 


*R.  D.  Leighty.  Optical  Power  Spectrum  Analysis  (OPSA)  (Report  No.  /,  Recording 
Optical  Spectrum  Analyser  System  Hardware) ,  U.  S.  Army  Engineer  Topographic 
Laboratories,  I-  o  rt  Belvoir,  VA,  Report  KTL-TR-74-1  1,  May  1975,  AD-A016  131. 

2R.  1).  Leighty,  and  (’■ .  K.  Lukes,  "‘Cloud  Screening  From  Aerial  Photography,”  Proc. 
dOtlt  Annual  Meeting  A  S  P  /  A  C  S  M  .  St.  L  <gu  is,  MO,  1  9  74  ,  p  p  309-332. 

1(i.  I- .  Lukes,  "Cloud  Screening  from  Aerial  Photography  Applying  Coherent  Optical 
Pattern  Recognition  Techniques,’’  Proc.  SPIF.,  Coherent  Optics  in  Mapping,  Vol.  45. 
I  974,  pp  265-272. 

4  (j .  k.  Lukes,  "Rapid  Screening  of  Aerial  Photography  by  OPS  Analysis."  Proc.  SPIF. 
O  ata  E  x  traction  prom  Film  ,  V  o  I.  1  1  7  ,  I  97  7  ,  p  p  8  9 -9  7 . 


14 


Spatial  Filter  Transform 
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(line  printer) 


C.  Software  processing  and  output 


FIGURE  4.  ROSA  Subsystems. 
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The  circular  sampling  aperture  a(x,  y)  of  radius  r  is  centered  on  the  optical  axis. 
When  illuminated  by  normally  incident  coherent  plane  waves  of  amplitude  B  and  wave¬ 
length  X,  the  aperture  has  a  transmission  function 


a (x. y) 


|B,  if  yj  x2  +  y2  <  r 


0,  otherwise 


(1) 


The  light  from  the  aperture  is  then  modulated  and  diffracted  by  the  transmission  func¬ 
tion  g(x,  y)  of  the  grating  transparency  to  yield  h(x,y).  such  that 


h  (x.  y)  =  a(x.y)  g(x.  y) 


C) 


The  diffraction-limited,  transform  lens  of  focal  length  f  and  diameter.  D, ,  where 
1)[  »  r.  then  images  the  Fraunhofer  diffraction  pattern  of  h(  x.y )  in  its  back  focal  plane. 
This  process  is  equivalent  to  a  two-dimensional  Fourier  transformation  F  (h  (x.  y)  ) 
except  fora  constant  phase  factor,  which  has  no  effect  upon  the  results.  Thus. 


F  (It  (x.  y  )> 


H  (u.v) 


1 

Xf 


// 

x.y 


a  (x.  y)  g  (x.  y) 


(31 


exp  |-j_ir(xu  +  yv)/(Xf)|  dxdy 

where  (u.  v)  represents  space  coordinates  in  the  detector  or  focal  plane.  To  express 
II  (u.v)  in  terms  of  spacial  frequencies  wu  and  wv.  use 

w„  =  and  wv  =  — (4) 


Titus,  H  (u /Xf.  v/XO  =  II  (wu,  wv). 

The  detector  is  responsive  only  to  light  intensity,  and  the  intensity  distribution 
at  the  photo  detector  is  equivalent  to  the  squared  modulus  of  the  light  amplitude,  or 

1  (u,  v)  =  j  II  (u.  v)  I-  (5) 

Now  consider  each  of  the  n  =  1 .64  photo  detector  elements  to  have  response 
K„  for  incident  light  energy.  The  signal  for  any  detector  element  is  then  a  function  of 
the  element  surface  area  and  the  responsivity  of  the  material.  At  this  point,  u  and  v 
can  be  expressed  in  polar  coordinates,  where  <>  =  s/  u-  +  v:  and  0  =  tan'1  (v/u). 
Therefore,  the  signal  from  any  detector  element  is 
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K 


n 


I  H  (p,  6)  I  2  p  dpd0,  n  =  1, 64 


(6) 


where  pn  and  0n  represent  both  the  lower  and  upper  limits  for  each  element  con¬ 
sidered  in  the  integral. 


Each  of  the  signals  next  passes  through  two  amplifier  stages,  the  first  has  constant 
gain  g^j  with  adjustable  bias  c“,  and  the  secoftd  has  variable  gain  g|j  with  an  ad¬ 
justable  bias  c|j.  The  net  amplification  beifip  • 


Sn  =  - 


4  b 
n 


(7) 


Signals  zn  at  the  output  of  these  ampli?  •«  .v  .  then 


i  A  (p,  0)  1 2  p  dpdfl 


(8) 


The  multiplexer  steps  through  each  of  the  n  channels  to  form  the  set  [zn], 
n  =  1 ,  64.  Next,  the  signals  pass  to  an  autoranging  amplifier  linked  to  a  binary-coded 
decimal  converter.  The  amplifier  has  four  gains  mp,  p  =  0,  3,  determined  by  the 
position  of  the  most  significant  digit  in  the  converter.  Final  output  of  the  amplifier  is 
then 


{  mP  z„  }  P  =  0,3  =  qn ,  n 


1, 64 


where 


m„  = 


/m„  =  2 

X 

10°, 

if 

104 

<z„  <  10s 

r  i 

II 

E 

X 

10', 

if 

I03 

<zn  <  104 

m2  =  2 

X 

102, 

if 

102 

<zn<103 

m,  =  2 

X 

103, 

if 

zn<102 

(9) 


(10) 


and  output  from  the  converter  is  then  the  three  most  significant  digits  and  the  decimal 
characteristic.  Thus, 

|Sn  }  =  (qn  /  10P,3-p),  p  =  0,3,  n  =  1,64  (11) 

The  encoder  then  reformats  {sn}  into  eight  groups  of  numbers  for  recording, 
along  with  a  six-digit  decimal  sample  descriptor  number  d  input  by  the  operator  from 
the  control  unit  and  card  number  i  =  1,8.  The  output  data  is  then  a  set  of  eight  card 
images 
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{s  }|  -  {d,  { Sj  =  ,  8  i}j  =  i  8  *  * 

from  which  the  sample  vector  S  =  {sn},  n  =  1,  64  is  available  for  the  computer 
pattern  recognition  efforts. 


Establishing  an  optimum  set  of  values  for  the  optical  and  electronic  parameters, 
i.e.  aperture  diameter,  transform  lens  focal  length,  amplifier  gains,  etc.,  to  obtain  the 
best  sample  vector  S  for  the  grating  patterns  to  be  investigated  is  dependent  on  know¬ 
ledge  of  the  effects  on  the  output  data  owing  to  variables  in  the  analytical  expressions 
just  discussed.  This  point  is  aptly  illustrated  when  reviewing  the  optical  power  spectrum 
analysis  literature  to  find  that  authors  do  not  recognize  the  effects  of  the  sampling 
aperture  on  the  output  data.  In  equation  (3),  the  sampling  aperture  a(x,y)  is  assumed 
to  have  a  spatially  uniform  or  constant  value  in  a  finite  area.  It  is  this  finite  area  that 
sets  the  integration  limits  with  the  result  that  equation  (3)  is  termed  the  superposition 
integral,  which  reflects  spatial  frequency  domain  convolution.  If  a(x,  y)  and  g(x,  y) 
have  respectively  the  Fourier  transforms  A  (u,  v),  and  G  (u,  v).  then  equation  (3) 
can  be  expressed  in  the  spatial  frequency  domain  by 

H(u,  v)  =  — -  A(u,  v)  **G(u,  v)  (13) 

Af 

where  the  double  astrisk  symbolizes  the  two-dimensional  convolution  operation,  which 
in  integral  form  is 


H(u,  v) 


A(a  — u,  P  -  v) 


G(u,  v) 


dad/3 


(14) 


Thus,  A(u,  v)  is  convolved  with  G(u,  v)  having  a  smoothing  effect,  and  the  amount 
of  smoothing  is  inversely  related  to  the  aperture  diameter. 


Aperture  effects  on  detection  of  gratings  are  of  considerable  importance  for  our 
problem.  Therefore,  we  will  provide  many  illustrations  related  to  diffraction  patterns 
resulting  from  circular  apertures.  These  illustrations  and  their  underlying  concepts  will 
be  referenced  frequently  in  the  following  discussions.  A(p)  is  represented  by  the  first- 
order  Bessel  function,  often  referred  to  as  the  Airy  function.  Hereafter,  this  will  be 
referred  to  as  the  Bessel  function.  Thus 


A(p) 


(rrpD/AQ 

(irpD/\f) 


(15) 
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whore  p  is  the  spatial  frequency  variable  and  1)  is  the  sampling  aperture  diameter. 
In  figure  6.  the  amplitude  distribution  for  circular  aperture  diffraction  patterns  is  shown 
for  two  values  for  aperture  diameter,  and  associated  with  each  of  these  patterns  are  three 
values  of  transform  lens  focal  length.  The  intensity  pattern  of  A(p), 


A  (p)> 


B27T-D4 


f  :j,  upD/xn  y 

L  OrpD/Xn  J 


(lb) 


is  pictured  in  figure  7  for  an  aperture  of  arbitary  diameter.  In  figure  8,  the  relation  of 
the  intensity  function  to  the  amplitude  function  without  weighting  is  indicated.  Note 
that  the  intensity  function  is  definitely  positive.  If  we  plot  the  intensity  function  semi- 
logarithmically  as  in  figure  lK  this  presentation  will  clearly  depict  the  intensity  function's 
relation  to  the  logarithmic  exposure  response  of  a  photographic  emulsion  as  shown  in 
figure  7.  If  the  aperture  diameter  is  decreased,  the  rings  in  figure  7  would  be  further 
spaced  as  shown  in  figure  6,  indicating  a  higher  degree  of  spatial  smoothing  in  equation 
(14).  This  spacing  is  significant  when  small  aperture  diameters  are  required  to  sample 
the  grating  maps  so  as  to  retain  the  boundary  definition  within  some  tolerable  minimum 
effective  ground-sampling  diameter.  The  scaling  effects  for  five  small  sampling  apertures 
and  two  different  transform  lens  focal  lengths  are  shown  in  table  I .  Note  the  decrease 
in  the  first  zero  crossing  (point  2)  with  increasing  aperture  diameter.  In  figure  10  is 
shown  the  fraction  of  the  total  energy  contained  within  circles  of  prescribed  radii  in  the 
diffraction  pattern  of  a  circular  aperture.  Approximately  84  percent  of  the  energy  lies 
within  the  central  lobe  of  the  Bessel  function.  Therefore,  to  minimize  smoothing  effects, 
this  first  zero  crossing  must  he  as  small  as  possible,  i.e.  the  diameter  of  the  sampling 
aperture  must  be  as  large  as  possible.  Also,  note  that  the  amplitude  values  for  the  tnax- 
imums  and  minimums  of  the  Bessel  functions  in  table  1  have  not  been  scaled  for  aperture 
area  according  to  equation  ( 1 5). 

f  rom  equation  (3)  the  effect  of  the  lens  focal  length  f  and  the  monochromatic 
wavelength  X  will  inversely  affect  the  amplitude  of  the  output  as  well  as  the  scale  of  the 
presentation  on  the  detector. 


The  optical  power  spectrum  has  polar  symmetry  as  indicated  in  equation  (5) 
since  I(u,  v)  =  I(-u,  -v).  This  redundancy  allows  use  of  the  ring  and  wedge  geometry 
in  design  of  the  ROSA  detector  without  any  lost  of  information. 


Detector  element  responses  Kn  are  a  function  of  the  respective  receptive  areas 
as  well  as  the  element's  responsivity.  Thus,  all  wedges,  being  of  the  same  area,  are  ex¬ 
pected  to  have  the  same  response  for  a  given  constant  responsivity;  whereas,  the  responses 
for  rings  will  increase  with  radius  from  the  center.  If  the  second  stage  amplifier  gain 
gjj  in  (7)  is  adjusted  so  that 


R I  R2  R3  R4  R5  R6 


ROSA  RINGS  (DETECTOR  MODEL  WRD-6420) 

FIGURE  6.  Amplitude  of  Circular  Aperture  Diffraction  Patterns. 
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TABLE  1.  Aperture-Diameter  Scaling  of 'Bessel  1  unction. 


LET  WHERE 

,  2J,(ffpD/X0  x  =  7rpl)/Xt' 

g(x)  =  — L -  x  x  r 

(jrpl)/Xl)  p  =  — :r-  ami  it  (  =  XJit 

IT  i)  ' 

T.  X1 

Then  p  =  (,  -p- 

DISTANCE  FROM  OPTICAL  AXIS,  p  in  nun. 


t  =  177. Smm.  X  =  0.6325  x  I  O''1  mm 

APERTURE  DIAMETER.  I)  in  mm. 


i 

C 

I 

Qlx) 

0.075 

0.10 

0.25 

0.50 

0.75 

1.00 

1 

0.0 

+  1.000 

0 

0 

0 

0 

0 

0 

■> 

1.220 

0 

1.830 

1.373 

0. 540 

0.275 

0. 1 83 

0.137 

3 

1.635 

-0.1323 

2.453 

1.870 

0.734 

0.368 

0.245 

0.187 

4 

2.233 

0 

3.350 

2.512 

1 .005 

0.502 

0.335 

0.251 

5 

2.670 

+0.0648 

4.010 

3.014 

1.206 

0.603 

0.402 

0.301 

6 

3.23S 

0 

4.858 

3.643 

1.457 

0.720 

0.486 

0.364 

7 

3.600 

-0.0400 

5.540 

4.1o2 

1 .665 

0.832 

0.555 

0.416 

S 

4.241 

0 

6.362 

4.772 

1 .000 

0.054 

0.636 

0.477 

t  =  1000.0 

mm 

X 

=  0.632 

8  x  10' 

mm 

APERTURE  DIAMETER,  D  in 

mm. 

i 

C  C?(  x ) 

0.075 

0.10 

0.25 

0.50 

0.75 

1.00 

l 

0.0  +1.000 

0 

0 

0 

0 

0 

0 

T 

1.220  0 

1 .020 

0.772 

0.300 

0.154 

0. 1 03 

0.077 

3 

1.635  -0.1323 

1.380 

1.035 

0.414 

0.207 

0.138 

0.104 

4 

2.233  0 

1.884 

1.413 

0.565 

0.283 

0. 1 88 

0.141 

5 

2.670  +0.0648 

2.260 

1 .605 

0.678 

0.330 

0.226 

0.170 

6 

3.238  0 

2.732 

2.040 

0.820 

0.410 

0.273 

0.205 

7 

3.600  -0.0400 

3.121 

2.341 

0.036 

0.468 

0.312 

0.234 

8 

4.241  0 

3.578 

2.683 

1 .073 

0.537 

0.358 

0.268 
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10.  t  raction  ot  the  Total  Inergy  C  ontained  Within  Circles  of  Pre¬ 
scribed  Radii  in  the  Diffraction  Pattern  of  a  Circular  Aperture. 
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then  the  responsivity  differences  of  detector  elements  is  ignored.  Although  it  gn  = 
constant,  i.e.  if  the  outputs  from  all  second-stage  amplifiers  are  equal,  then  detector 
element  areas  are  normalized.  However,  the  flexibility  exists  with  software  to  adjust 
data  for  desired  area  response,  thus,  the  former  approach  is  generally  used  m  practice. 


B.  Gratings  in  OPS  Systems  •  A  transmission  grating  will  be  defined  as 
a  one-dimensional  repeating  optical  transmission  function.  For  example. 

g  ( x ,  y )  =  g  ( x )  =  A  cos(2rrx/p)  (lb) 

is  a  cosine  function  of  amplitude  A  in  the  x-direction  that  repeats  every  p  units  of 
x,  and  the  value  of  the  function  at  any  point  x  is  the  same  for  all  values  of  y.  The 
spatial  frequency  of  a  grating  in  line  pairs  per  millimeter  (mm)  is  the  reciprocal  of  the 
grating  period  in  millimeters,  i.e. 

w  =  I  /  p  (Id) 

The  Fourier  transform  of  g  (x)  is  G  (w),  which  will  be  symbolized  as  shown  in  equation 
(3),  i.e.  F  Ig  (x)|  =  G  (w).  As  in  equation  (18)  where  the  grating  can  be  expressed 
in  one  dimension,  it  is  frequently  advantageous  to  present  expressions  as  F  |g  (x.y)  = 
g  (x)|  =  G  (u,  o)  G  (p,  0).  where  p  is  a  radial  distance  and  0  is  a  direction  angle. 

Since  the  grating  transform  for  a  specific  grating  will  have  0  =  a  constant,  we  can  then 
talk  in  general  of  G  (p,  0 j)  =  G  (p).  By  equation  (4),  we  have  w  =  p/Xf,  thus 
G  (p/AO  =  G  (w). 

As  shown  in  figure  1 1,  the  Fourier  transform  for  a  cosine  grating  is  equal  to  two 
delta  functions,  5(w  +  I/p),  with  amplitudes  of  A/2  positioned  at  +l/p  in  the  spatial 
frequency  presentation.  This  results  from  the  equality 

cjw<>*  +  e-jw0x 

COS(W„  x)  =  - - - ,j2  =  -1 


and  since  F(eJW°']  =  5  (w-w0),  where  wQ  =  1/p  in  lp/mm.  we  have  for  equa¬ 
tion  (18) 


A  few  words  about  delta  functions  is  necessary.  A  delta  function  6  (t),  also 
known  as  a  unit  impulse  function,  may  be  defined  in  many  ways.  Here  the  following  definition 
will  be  used: 


10,  if  t  =  0 

“,  if  t  =  0 


i 


6  (t)  dt 


e 


i 


6  (t)  dt  =  1,  e  >  0 


(21) 


The  first  relation  indicates  that  5  (t)  is  zero  except  at  t  =  0,  where  it  becomes  infinite 
in  such  a  way  that  the  second  relation  is  satisfied.  Delta  functions  have  many  useful 
mathematical  properties,  which  we  will  only  reveal  as  needed  to  clarify  the  concepts 
under  discussion. 


Delta  functions  in  optics  are  point  sources  of  light;  and  since  the  cosine  grating 
has  only  two  delta  functions,  each  of  which  has  infinite  resolution  in  the  frequency 
domain,  it  would  appear  that  this  grating  would  be  ideal  for  OPSA  applications  such 
as  we  are  investigating.  Indeed,  if  we  had  means  for  fabricating  cosine  gratings  of  selected 
spatial  frequencies  with  the  desired  spatial  boundary  conditions,  it  would  be  ideal; 
unfortunately  we  do  not. 

A  square  wave  grating  would  be  much  simpler  to  realize,  except  that  we  cannot 
obtain  negative  light  values  with  optical  transparencies.  We  can  however  use  a  'Ronchi' 
grating,  which  is  a  square  wave  grating  with  a  positive  A/2  bias. 
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Another  unique  characteristic  of  a  Ronchi  grating  is  that  the  grating  width,  d, 
is  one-half  of  the  grating  period,  p.  The  disadvantage  of  using  Ronchi  gratings  instead 
of  cosine  gratings  is  that  harmonic  delta  functions  are  found  in  the  Fourier  transform, 
as  shown  in  figure  11.  These  harmonics  decrease  in  amplitude  with  increasing  spatial 
frequency  and  alternate  plus  and  minus  in  amplitude.  To  obtain  the  functional  expression 
for  Ronchi  grating  transforms,  we  use  the  Fourier  series  representation  for  a  square 
wuvetrain  with  a  constant  bias.  Using  the  Fourier  series  results  in  an  infinite  number  of 
coefficients.  An  infinite  number  of  cosine  waves  are  required  to  realize  perfect  square 
waves,  but  if  we  select  the  first  n  coefficients,  where  n  is  a  modest  number,  we  can 
transform  each  Fourier  series  term  and  the  result  is  that  each  Fourier  series  term  has 
its  unique  delta  function  in  the  spatial  frequency  domain. 


By  using  Ronchi  gratings,  we  trade  simplicity  of  grating  fabrication  for  complexity 
in  detection  of  the  grating  transforms.  We  will  show  that  the  detection  of  ideal  Ronchi 
gratings  can  be  made  at  sufficiently  high  signal-to-noise  ratios  to  yield  unique  signals  for 
OPSA  applications. 


('.  (i rating  and  Aperture  Transform  Convolution.  •  liquation  (14)  is  the 
convolution  integral  for  a  general  signal  function.  Cl  (u,  v).  and  a  general  aperture 
function,  A  (u.  v).  In  this  section  we  will  be  more  specific  and  use  grating  functions 
and  small  circular  sampling  aperture  functions.  The  objective  is  to  indicate  the  nature 
of  the  resulting  function  indicated  by  equation  (14)  when  parameter  values  of  the 
aperture  and  grating  function  are  varied. 


To  introduce  the  concepts  with  minimum  complexity,  we  will  first  assume  that 
a  cosine  grating  is  illuminated  by  a  collimated,  coherent  light  beam  passing  through  a 
circular  sampling  aperture.  We  will  rewrite  equation  (14)  in  one-dimensional  form  as 


II  (w) 


A  (a)  (I  ( w  -  a) 


da 


Now  substitute  equations  (15)  and  (  20) 


II  (w) 


1  r"  -J,  (irnD/Xf)  A 

77  I  : - ' - T777 - 7  I  6(w-  l/p-a)  +  S(w+  l/p-a)|  da 

At  J  4  (7ra()/Al)  2 
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With  the  relation. 


Then, 


F  (w)  *  6(w-wo)  =  J  (a-  w^)  F  (a)  da  =  F(wo) 


1  ABttD2  f  2J . (7r(w-l /p)D/Xf)  2J,(jr(w+l/p)DA0  | 

H  (w)  =  —  -  — ! -  +  — -  I  (22 

Xt  8  L  (ir(w-l/p)D/XO  (rr(w+l/p)D/XO  J 


2) 


Hxpression  (22)  tells  us  that  convolving  an  optical  cosine  grating  transform  with 
the  first-order  Bessel  function,  which  is  the  transform  of  the  aperture  function,  is  equiva¬ 
lent  to  creating  two  new  first-order  Bessel  functions  with  one  positioned  at  the  spatial 
frequency  1/p,  and  the  other  at  -1/p.  The  Bessel  functions  will  have  their  same  spatial 
frequency  distribution  as  that  of  the  aperture  function,  but  will  be  scaled  by  a  factor 
of  one-half  in  amplitude.  Another  important  revelation  of  expression  (22)  is  the  poten¬ 
tial  for  significant  interaction  between  the  two  Bessel  functions  resulting  from  the  con¬ 
volution.  We  have  seen  from  figure  6  and  table  1  how  the  Bessel  function  is  scaled  in 
spatial  frequency  inversely  with  the  diameter  of  the  sampling  aperture.  Then,  for  very 
small  sampling  apertures  and  low  spatial  frequency  gratings,  one  could  expect  to  have 
overlapping  Bessel  functions  that  add  algebraically  to  yield  a  net  amplitude  spectrum 
somewhat  different  than  that  of  the  individual  Bessel  functions.  We  will  show  examples 
of  these  implications  in  the  simulations  to  follow. 


Using  the  results  of  expression  (22)  as  a  guide,  we  can  easily  extrapolate  to  a 
simple  expression  for  an  ideal  Ronchi  grating  illuminated  by  collimated,  coherent  light 
passing  through  a  circular  sampling  aperture.  The  major  differences  between  the  ex¬ 
pressions  that  result  from  the  differences  between  the  cosine  transform  function  and  the 
Ronchi  grating  transform  function  are  as  indicated  earlier,  and  as  shown  in  figure  11. 
The  infinite  expression  for  Ronchi  grating  convolution  is  then  more  complicated  in 
appearance  and  is  given  by 
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This  convolution  is  indicated  in  figure  1 2.  Note  the  net  amplitude  spectrum 
resulting  front  the  individual  overlapping  Bessel  amplitude  spectra.  Since  only  the  in¬ 
tensity  spectrum  is  detectable,  the  bottom  of  figure  I  2  shows  the  intensity  spectrum 
to  be  the  magnitude  squared  of  the  net  amplitude  spectrum,  as  indicated  in  equation  (6). 
Thus, 


I  (w)  =  IH  (w)  P 


(24) 


To  understand  the  intensity  function  of  equation  (24)  in  terms  of  the  individual 
Bessel  function  in  equation  (23),  we  will  generalize  (23)  to 

H(w)  =  K  [  A,  J  (w)  +  A,J(w  +  w,)  +  A2J  (w  -  w, ) 

-A3J(w  +  3w,)  -  A3J(w-3w])  +  ...|  (25) 


The  first  term  within  the  parenthesis  is  a  Bessel  function  centered  on  the  optical 
axis,  the  second  and  third  terms  represent  the  positive  and  negative  first  harmonics, 
respectively,  and  the  fourth  and  fifth  terms  represent  the  positive  and  negative  third 
harmonics.  Higher  order  terms  have  been  omitted  from  equation  (25). 

When  ( -5)  is  squared ,  as  in  ( 24),  interaction  terms  result. 

H2  (w)  =  K2  [A2J2  (w)  +  A2J2  (w  +  w()  +  A2J2  (w  -  w( ) 

-  A2  J2  (w  +  3w[ )  -  A2J2  (w  +  3w  )  +  .... 

+  A,  A2J(w)J(w  +  w|)  +  A|A2J(w)J(w-w  ) 

-  A|A3J(w)J(w  +  3w))  -  A(  AjJ  (w)  J  (w  -  3W|  ) 

-  A2A3J(w  +  w1)J(w+3w1)  -  A2A3J(w-w1)J(w-3w1) 
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FIGURE  12.  Convolving  Amplitude  Spectra  From  the  Sampling  Aperture 
and  Grating  to  Yield  a  Detectable  Intensity  Spectrum. 


t 


The  significance  of  the  interaction  terms  is  dependent  upon  the  scaling  of  the 
Bessel  function,  as  dictated  by  the  size  of  the  sampling  aperture  and  the  focal  length 
of  the  transform  lens,  and  upon  the  positioning  of  the  functions  as  determined  by  the 
Roik'  i  grating  frequency,  w( .  If  the  sampling  aperture  is  large,  the  interaction  terms  will 
be  insignificant;  but  when  the  sampling  aperture  is  small  and  the  grating  frequency  is 
small,  then  some  of  the  interaction  terms  can  be  significant.  This  will  be  discussed  later. 


It  should  be  noted  that  (23)  implies  that  the  Ronchi  grating  is  a  real  function, 
i.e.,  has  no  phase  components.  In  practice  this  is  never  fully  realized,  and  the  effects 
include  even-order  harmonics.  This  will  be  discussed  subsequently. 


D.  Geometry  of  ROSA  Detectors.  •  Equation  (23)  indicates  the  nature  of 
the  amplitude  signal  on  the  detector.  Here  we  describe  the  geometrical  nature  of  ROSA 
detectors  that  will  partition  or  sample  the  incident  grating  spectra,  and  the  following 
section  will  discuss  the  sampling  of  grating  spectra. 

Two  ROSA  detector  models  are  available  from  Recognition  Systems  Inc.,  each 
has  32  ring  elements  and  32  wedge  elements;  however,  the  detectors  have  some  major 
differences.  The  detection  of  Ronchi  gratings  was  simulated  for  both  detectors  with 
slightly  different  results  that  are  better  understood  with  a  knowledge  of  the  detectors 
and  their  differences. 


The  older  detector.  Model  WRD-6420,  lias  two  20°  scratch  wedges  (usually 
oriented  normal  to  the  direction  of  roll-film  travel  in  a  manner  such  that  emulsion 
scratches  will  not  yield  false  diffraction  signals),  and  the  newer  detector.  Model  WRD- 
6400.  has  no  scratch  wedges.  The  face  plates  of  these  detectors  is  shown  in  figure  13. 
Table  2  lists  the  dimensions  of  the  ROSA  Model  WRD-6420  detector  and  table  3  presents 
similar  data  for  the  Model  WRD-6400  detector.  In  both  tables,  the  data  is  presented  for 
each  detector  ring  element.  The  least  and  maximum  radius  and  their  radial  difference  give 
the  dimensions  of  the  concentric  detector  rings  that  integrate  incident  energy.  For  each 
ring,  the  mid-radius  is  a  key  dimension  in  the  later  simulations  because  the  gratings  will 
be  designed  to  have  their  first  harmonic  positioned  at  the  mid-radius.  The  tables  indicate 
the  spatial  frequencies  associated  with  the  mid-radii  for  four  different  focal  length  trans¬ 
form  lenses.  In  figure  14,  the  relation  is  plotted  between  radial  distance  and  spatial 
frequency  in  the  back  focal  plane  of  these  four  lenses  as  given  by  equation  (4).  The 
remaining  dimension  shown  in  tables  2  and  3  concerns  the  wedge  width  at  niid-radius. 
This  dimension  will  be  used  to  analyze  grating  direction  detection  in  the  simulations 
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TABLE  2.  KOSA  Detector.  Model  WRD-6420 
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FIGURE  14.  Spatial  Frequency-Radial  Distance  Relation  in  the  Back 

Focal  Blanc  of  a  Transform  Lens. 


The  wedge  angles  lor  t he  WRl)-b420  ileleelor  are  5.000°,  and  the  wedge  angles 
tor  the  WRI)-b400  deteetor  are  5.b25J:  therefore .  the  wedge  widths  will  be  different 
for  each  deteetor. 


Note  that  the  inner  radius  for  WRl)-b420  detector  wedges  is  1.17b  millimeter. 
Therefore,  only  ring-wedge  combinations  for  rings  7  through  52  are  available  for  grating 
detection,  for  the  Model  WRl)-b400  detector,  wedges  have  an  inner  radius  of  0.544 
millimeter,  except  for  wedges  lb  and  17.  which  have  an  inner  radius  of  0.660  millimeter 
corresponding  to  rings  b  and  7  respectively.  The  outer  radius  of  the  wedges  for  Model 
WRl)-b420  detectors  is  1 5.3b7  millimeters:  whereas,  the  maximum  radius  for  Ring  32 
is  15.875  millimeters.  This  will  have  significance  in  later  simulations,  for  the  WRI)-b400 
detector,  the  outer  wedge  radius  and  the  maximum  radius  for  Ring  32  are  all  equal  to 
1 5.000  millimeters. 


In  figures  15  and  1  (>.  linear  and  semi-log  plots  of  the  mid-radii  for  the  two  de¬ 
tectors  are  presented.  In  general,  figure  15  indicates  that  the  WRl)-b400  detector  has 
greater  frequency  domain  resolution  than  the  WRl)-b420  detector  up  to  about  Ring  20 
and  less  resolution  at  higher  spatial  frequencies,  figure  lb  indicates  that  the  rings  of 
detector  WRl)-b400  are  spaced  logrithmieully  beyond  Ring  b.  In  figures  17  and  18, 
linear  and  semi-log  plots  of  the  ring  widths  for  the  two  detectors  are  presented.  Here 
figure  17  shows  the  ring  widths  for  the  WRI)-b420  detector  to  increase  in  a  linear  man¬ 
ner:  whereas,  those  of  the  WRl)-b400  detector  are  shown  to  he  nonlinear.  In  the  semi-log 
plot  of  figure  18.  both  are  quasi  semi-log  in  nature,  but  with  different  slopes. 
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l  'CjlIRt  15.  Linear  Plot  of  Mid-Ring  Radii  for  Two  ROSA  Detectors. 
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FKiliRK  17.  l  inear  Plot  of  Detector  Ring  Width  lor  Two  ROSA  Detectors. 
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F.  ROSA  Sampling  of  Grating  Spectra.  •  In  the  ROSA  characterisation 
section,  we  analytically  outlined  the  general  optical  and  electrical  signal  tlow  through 
the  ROSA  system.  The  next  section  introduced  gratings  in  OPS  systems,  and  the  fol¬ 
lowing  section  dealt  with  grating  and  aperture  convolution.  Then  the  previous  section 
introduced  the  geometry  of  the  ROSA  detectors.  In  this  section,  we  integrate  the  pre¬ 
vious  material  and  discuss  Ronchi  grating  sampling  with  ROSA  detector  ring  and  wedge 
elements. 


For  this  discussion,  we  assume  the  condition  of  equation  (17)  to  be  valid,  i.e., 
the  gain  and  bias  of  the  amplifiers  are  adjusted  to  normalize  responsivity  differences  of 
the  detector  elements  and  also  normalize  detector  element  areas.  This  condition  implies 
a  point  source  of  light  imaged  on  any  detector  element  will  yield  the  same  output.  This 
simplifying  assumption  in  no  way  restricts  the  generality  of  the  approach,  as  a  review  of 
the  ROSA  characterization  section  will  indicate  that  any  set  of  amplifier  gains  and  biases 
can  be  used  to  compute  any  other  set,  if  so  desired. 


Now,  consider  that  each  ring  element  of  the  detector  has  a  finite  width  with  an 
inner  radius  and  outer  radius  and  that  each  ring  element  is  separated  by  a  nonsensitive 
band.  We  can  use  these  radii  for  each  ring  for  the  limits  of  integration  in  equation  (6) 
for  determining  the  signal  from  any  intensity  spectrum.  This  can  be  visualized  by  con¬ 
sidering  the  one-dimensional  intensity  spectra  at  the  bottom  of  figure  12  with  the  de¬ 
tector  ring  radii  converted  to  their  spatial  frequency  equivalents  for  a  given  scaling 
factor  of  the  transform  lens  focal  length  and  wavelength  of  light.  The  signal  for  each 
detector  element  is  then  proportional  to  the  area  of  the  curve  between  the  ring's  inner 
and  outer  radii. 


We  can  see  from  figure  I  that  there  may  be  one  or  more  Bessel  functions  within 
the  area  of  a  wedge,  depending  upon  the  spatial  frequency  of  the  Ronchi  grating  illu¬ 
minated  by  the  aperture.  Suppose  one  side  of  the  intensity  spectrum  in  figure  12  re¬ 
presents  the  center  line  of  a  wedge;  the  signal  for  the  detector  wedge  element  can  be 
visualized  by  considering  the  area  under  the  curve.  Depending  upon  the  diameter  of  the 
sampling  aperture  and  scaling  by  the  transform  lens  focal  length,  the  adjacent  detector 
wedge  elements  will  have  detectable  signal  levels  owing  to  the  portions  of  the  Bessel 
functions  imaged  onto  their  surfaces. 


If  we  can  uniquely  fabricate  ideal  Ronchi  rulings  at  2b  selected  grating  periods 
such  that  with  a  selected  focal  length  of  a  transform  lens  a  unique  signal  will  result  for 
each  ring  element  of  the  detector  between  7  and  32.  and  if  we  further  have  the  capability 
of  fabricating  these  Ronchi  gratings  with  5-degree  differences  in  direction,  we  then  have 
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the  potential  capability  of  encoding  832  separate  types  of  information  for  the  WRD-6420 
detector.  In  the  following  section,  we  will  simulate  the  grating  detection  process  for 
both  detectors  to  show  that  unique  signals  can  be  detected,  even  for  small  sampling 
apertures. 


A  few  comments  on  the  sampling  at  grating  boundary  locations  in  OPS  systems 
are  warranted  in  discussing  ROSA  grating  sampling.  Let  us  begin  with  equation  (3), 
the  basic  expression  for  the  grating  illuminated  by  coherent  light  passing  through  a 
sampling  aperture.  Suppose  the  grating  map  g  (x,  y)  is  such  that  g(x,  y)  =  gj  (x,y)  + 
g,  (x,  y),  where  g,  (x,  y)  and  g,(x.  y)  are  nonoverlapping,  contiguous  grating  func¬ 
tions.  For  the  moment,  suppose  the  sampling  aperture  diameter  is  very  large:  therefore, 
the  aperture  convolution  effects  can  be  neglected,  e.g.  if  D  =  50.0  mm  then  the  central 
lobe  of  the  Bessel  function  has  a  diameter  of  0.0055  mm.  as  indicated  in  table  1 
for  a  transform  lens  focal  length  of  177.8  mm  and  a  wavelength  of  light  of  0.6328  x 
1 0'-1  mm. 


If  the  illumination  is  uniform  over  the  aperture  and  the  aperture  is  located  im¬ 
mediately  adjacent  to  the  transparency  along  the  optical  axis,  then 

g  (x.  y)  =  c,  g,  (x.  y)  +  c,  g2  (x,  y)  (27) 

where  c,  and  c,  represent  the  percentage  of  the  aperture  illuminated  by  g,  (x.  y) 
and  g,  (x,  y),  respectively.  If  the  respective  Fourier  transforms  of  g,  (x,  y)  and 
g,  (x,  y)  are  C»,  (u,  v)  and  G,  (u,  v),  then  the  linearity  property  of  the  Fourier  trans¬ 
formation  indicates 


[Cj  g,  (x,  y)  +  c2  g2  (x,  y)|  =  — -7  (c,  G,  (p.  0, )  +  c2  G,  (p,  0,)) )  (28) 

A! 

when  expressed  in  polar  coordinates  for  easier  visualization.  In  this  expression,  the  delta 
functions  of  G(  (p.  0{)  and  G,  (p,  0,)  are  separately  collinear  along  their  respective 
directions  of  0,  and  0,.  The  intensity  function  of  equation  (28),  however,  indicates 
the  possible  existence  of  interaction  between  G,  (p,  0()  and  G,  (p,  0,),  as  shown  by 

l(p,0)  =  — i-  |  c,  (>,  (p,  0( )  +  c,  G,  (p,  0,)  1  2 

=  TT77  (cfG^p.O,)  +  2c,  G  (p,  0  )  c,  G,  (p,  0  ) 

A"  I- 

+  c22C,](p.02))  (29) 
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When  0,  differs  from  0,,  the  interaction  term  concerns  only  the  della  functions  at 
the  origin.  When  0,  =  0\ ,  there  will  be  no  additional  interaction  unless  the  periods 
of  (x,  y)  and  g->  (x,  y)  are  integer  multiples,  which  as  table  2  indicates  is  not  likely 
with  the  ROSA  Model  WRD-6420  detector  when  the  Ronchi  gratings  are  properly  de¬ 
signed  . 


Now,  let  us  generalize  the  discussion  for  any  sampling  aperture  diameter  so  as 
to  consider  convolution  smoothing  and  interaction'in  the  spatial  frequency  domain. 
Expression  (28)  is  now  rewritten  as 


ja(x.y)  |c,  g,  (x,y)  +  c2g2(x,y)j^  =  -  |a(p,0)  **  ^c,G,(p,01)  + 

c2G2(p.0,))]  =  jL  (a (p,0)  **  c,  G,  (P.0,)  +  A  (p.0)  ** 

c,  G,  (p.02)) 


(30) 


The  intensity  function  for  equation  (30)  is 


'A(p,0)  **  Cj  G  (p,  0  j )  2  +  2  |a  (p,  0)  **  c,  G,  (p,0,) 

(31) 

A  (p,  0)  **  c2  G2  (p.02))+  (  A(p,0)  **  c2  G2  (p.02)j2 

The  first  and  third  terms  of  equation  (31)  are  additive  where  they  overlap  in  (p,  0) 
space  owing  to  the  convolution  broadening.  The  second  term  produces  a  multiplicative 
function  at  the  origin  that  is  added  to  those  of  the  first  and  third  terms.  Additionally, 
the  multiplicative  interaction  of  the  amplitude  functions  in  (p,  0)  space  occurs  where 
the  functions  of  terms  one  and  three  overlap.  Although  this  interaction  could  be  con¬ 
siderable  for  very  small  sampling  apertures,  in  the  Systems  Concepts  Section  a  number 
of  strategies  will  be  discussed  for  uniquely  determining  grating  frequency  and  direction 
when  the  sampling  aperture  illuminates  more  than  one  grating  i unction. 

This  discussion,  although  being  somewhat  ^eneral  in  interpretation  to  the  real 
world,  pertains  mainly  to  the  idealized  or  theoretical  Ronchi  spectra.  Real  world  spectra 
will  contain  even-order  harmonics  with  amplitudes  somewhat  proportional  to  their  de¬ 
viation  from  the  theoretical.  The  real  world  effect  on  the  theoretical  data  will  result  in 
the  potential  for  overlapping  spectra  of  odd  and  even  harmonics  with  an  associated  loss 
in  spatial  frequency  discrimination. 
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IV.  SIMULATING  OPS  DETECTION 
OF  RONCHI  GRATINGS 


A  FORTRAN  IV  program,  WRDET,  has  been  prepared  to  simulate  the  concepts 
of  the  preceding  section,  wherein  ideal  Ronchi  gratings  of  selected  spatial  frequencies 
are  illuminated  with  collimated  light  passing  through  a  circular  sampling  aperture.  The 
intensity  spectra  derived  from  the  convolved  Fourier  transforms  for  these  two  functions, 
i.e.  the  grating  and  the  aperture  functions,  as  given  in  equation  (26),  is  then  partitioned 
and  summed  for  the  appropriate  ring  and  wedge  elements  of  the  ROSA  detector  con¬ 
figuration.  The  critical  variables  tested  were  sampling  aperture  diameter,  Ronchi  grating 
spatial  frequency,  and  focal  length  of  the  transform  lens. 


In  this  Chapter,  we  will  describe  the  major  features  of  WRDET  and  present  simu¬ 
lation  results. 


A.  Program  WRDET.  •  Program  WRDET  simulates  detection  of  ideal  Ronchi 
optical  transmission  gratings  in  an  optical  power  spectrum  analyzer  using  a  ROSA  de¬ 
tector.  For  a  given  transform  lens  focal  length  and  diameter  of  the  sampling  aperture, 
the  OPS  intensity  distributions  for  RONCHI  gratings  having  spatial  frequencies  corres¬ 
ponding  to  mid-ring  radii  of  the  ROSA  ring  detector  elements  were  partitioned  according 
to  the  ROSA  wedge  and  ring  detector  elements  and  integrated  for  selected  wedge  and 
ring  detector  elements.  All  calculations  were  accomplished  for  the  HeNe  laser  light  wave¬ 
length  of  0.6328  micrometer  (pm). 


The  discussion  will  follow  the  How  chart  of  the  major  parts  of  Program  WRDET 
given  in  figure  19.  A  listing  of  the  program  is  included  as  appendix  A.  Initialization  of 
Program  WRDET  involves  setting  the  values  for  the  sampling  aperture  diameter,  the 
transform  lens  focal  length,  the  Ronchi  grating  spatial  frequencies  corresponding  to  the 
mid-ring  radii  for  the  selected  ring  detector  elements,  the  inner  and  outer  radii  for  the 
ring  and  wedge  detector  elements,  and  a  list  of  the  first  tour  Bessel  zero  coefficients. 
An  optical  scaling  factor  for  the  Bessel  function  is  computed,  and  this  factor  is  used  to 
compute  the  radii  of  the  Bessel  rings  in  the  spatial  frequency  plane,  as  indicated  in 
table  1 . 


The  main  DO  loop  of  Program  WRDET  then  increments  through  each  selected 
grating  frequency.  For  each  selected  spatial  frequency,  the  amplitude  and  radial  distances 
are  computed  for  each  delta  function  falling  on  the  detector,  as  indicated  in  figures 
I  I  and  I  2. 
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For  each  mid-ring  spatial  frequency,  the  integrated  OPS  is  computed  for  the 
selected  ring  Rj  and  its  adjacent  rings  RJ+1  and  Rj_j  .  Similar  computations  are 
made  for  a  wedge  Wj  and  its  neighbor  WJ+1  when  the  grating  is  oriented  so  that  the 
Ronchi  grating  deltas  are  collinear  along  the  bisector  of  the  wedge  Wj.  These  com¬ 
putations  are  made  in  Subroutine  A1RYD,  which  is  outlined  in  the  flowchart  in  Figure  20. 
This  subroutine  integrates  the  two-dimensional  diffracted  energy  incident  on  a  wedge 
or  ring  detector  element  in  the  spatial  frequency  plane  by  Simpson’s  Rule  using  polar 
coordinates.  Thus,  the  subroutine  performs  the  integration  of  equation  (26)  where  the 
definite  integrals  have  limits  set  by  the  indicated  wedge  or  ring  element. 


When  Subroutine  A1RYD  is  called,  the  radial  and  angular  sampling  parameters 
are  computed.  Next,  a  table  of  normalized  first-order  Bessel  amplitudes  are  computed 
for  the  appropriate  scaling  parameter,  as  indicated  in  figure  8.  Then,  nested  integration 
DO  loops  increment  through  all  angle  and  radial-distance  sample  points.  For  each  sample 
point,  the  radial  distance  between  that  point  and  the  position  of  each  Ronchi  delta 
function  on  the  detector  is  used  to  determine  the  scaled  Bessel  amplitudes.  The  sum  of 
these  amplitudes  for  each  point  are  then  squared  to  obtain  the  intensity  estimate. 
Simpson's  integration  rule  is  applied  for  the  radial  and  angular  dimensions  to  provide 
a  number  representing  the  detector  element  signal. 


The  major  assumption  associated  with  Program  WRDFT  involves  the  maximum 
radius  of  the  first-order  Bessel  function  used  in  the  computation  of  diffraction  ampli¬ 
tudes.  The  first-order  Bessel  function  is  infinite  in  extent  and  to  make  computation 
tractable,  yet  meaningful,  required  truncation  of  the  Bessel  function.  As  indicated  in 
figure  10,  the  function  has  most  of  its  energy  within  the  central  lobe.  83.8  percent,  and 
the  total  energy  contained  by  the  Bessel  function  radii  within  the  second,  third,  and 
fourth  zeros  are  91 .0.  93.8.  and  95.3  percent,  respectively.  In  table  1 .  the  radial  distances 
for  these  zero  crossings  are  shown  for  six  aperture  diameters  and  two  transform  lens 
focal  length  values.  A  realistic  and  effective  compromise  for  the  Bessel  function  trun¬ 
cation  radius  would  lie  the  fourth  zero  radius. 


The  line-printer  output  from  Program  WRDF.T  needs  to  he  explained.  The  results 
of  the  simulation  data  are  presented  in  appendix  B  and  will  be  discussed  in  detail  sub¬ 
sequently.  Kach  simulation  run  is  on  a  seperate  sheet  having  a  header  indicating  pro¬ 
gramed  values  for  transform  lens  focal  length,  sampling  aperture  diameter,  radius  to  the 
fourth  Bessel  zero,  optical  wavelength,  RSI  detector  model,  and  number  of  sampling 
divisions  in  the  Bessel  function  to  the  fourth  zero.  The  body  of  the  printout  lists  the 
integrated  values  tor  equation  (26)  lor  detector  rings  R,  t ,  R(,  and  R|+)  and  detector 
wedges  W(  and  W|+  (  for  each  mid-ring  grating  spatial  frequency . 
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Loop  for  Each 
0  Sample 


Loop  for  Each 
R  Sample 


Loop  for  Each 
Ronchi  Delta 


FIGURE  20.  Flow  ('hart  of  Subroutine  AIR  YD. 


It  is  important  to  note  that  the  simulation  values  have  not  been  scaled  hv  the 


factor 

>s,  =  K- 

as  shown  in  equation  (23)  and  equation  (26).  In  t  >le  4.  the  values  for  the  scaling  factor 
are  indicated  for  four  transform  lens  focal  lens  i  ,i..J  four  aperture  diameters  when 
B  =  1.0.  To  adjust  the  light  input  amplitude  B  so  as  to  maintain  normalized  output 
from  the  values  in  table  4,  we  have 


JM2: 

4Xf 


4Af 

trl)' 


(33) 


and  the  expected  range  of  B.  for  the  simulations  to  be  presented,  is  well  within  oper¬ 
ational  optical  capabilities. 


B.  Simulation  Results.  •  Program  WRDIT  was  used  to  provide  simulation 
data  for  the  OPS  detection  of  Konchi  gratings  with  a  ROSA  detector.  The  simulations 
are  grouped  into  three  data  sets  and  summarized  with  their  simulation  conditions  in 
table  5.  Data  Set  1A  simulates  detection  conditions  for  the  WRD-6420  detector  with 
different  transform  lens  focal  lengths  and  sampling  aperture  diameters,  and  Data  Set  IB 
simulates  the  WRD-6400  detector  with  similar  conditions.  Data  Set  2  simulates  data 
from  uniform  Ronchi  grating  deformations  for  the  WRD-6420  detector,  i.e.  the  Ronchi 
grating  frequency  is  displaced  some  percentage  from  the  mid-ring  grating  frequency  in 
Data  Set  2A  and  the  Ronchi  grating  direction  is  displaced  by  an  angular  factor  from  the 
mid-wedge  angle  in  Data  Set  2B.  Data  Set  3  simulates  a  non-existant  detector  with  equal 
ring  widths  (Data  Set  3A)  and  with  modifications  of  equal  ring  width  (Data  Set  3B). 
Table  5  lists  the  illustration  number  for  the  graphical  presentation  of  the  data  and  the 
appendix  B  page  of  the  data  listings. 
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TABLE  4.  Intensity  Scaling  Eactors. 
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where  B=1.0  and  x=0. 6328x10” 3mm 


Aperture 


Diameter 

f  =  100.0  mm 

f  =  177.8  mm 

f  =  622.3  mm 

f  =  1016 

0.1 

0.01540 

0.99487 

0.00040 

0.00015 

0.25 

0.60174 

0.19035 

0.01554 

0.00583 

0.5 

9.62780 

3.04553 

0.24861 

0.09327 

1.0 

154.4472 

48.72853 

3.97784 

1.49231 
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TABLfc  5.  Summary  of  Simulation  renditions  with  Illustration 

and  Data  Page  References. 

(all  measurements  in  mm) 
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TABLE  5.  Continued. 


Data  Set  2B: 

WRD-6420  Detector,  D  =  0.1  mm,  f  =  100.0  mm,  Grating  angle  =  wedge  center 

+  factor  angle 
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Data  Set  3B: 

Modified  equal  ring  width  detector,  D  =  0.1  mm,  f  =  100.0  mm 
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1.  Data  Set  1A  -  Detector  WRD-6420  Simulations.  Figures  21  through 
27  contain  plots  of  the  simulation  data  for  Data  Set  1A.  The  purpose  of  this  data  set 
is  to  show  the  nature  of  the  detector  element  responses  and  S/N  for  the  WRD-6420 
ROSA  detector  with  different  transform  lens  focal  lengths  ( 1 77.8  and  100.0  mm)  and 
sampling  aperture  diameters  (0.5.  0.25.  0.10.  and  0.075  mm).  In  reviewing  this  data 
set,  we  will  first  look  at  the  gross  comparison  of  the  curves  and  then  at  the  details  of 
curve  configuration  and  curve  differences.  The  WRD-6420  dimensional  data  in  table  2 
and  figures  2  and  18  will  be  referred  to  frequently  in  the  detailed  discussions. 


For  a  transform  lens  focal  length  of  177.8  mm.  we  note  that  the  S/N  for 
the  rings  and  wedges  (R|/Rt.1.  or  Rj/Rj+i.  and  Wj/WJ+1)  increases  sharply  for  the 
0.5  mm  sampling  aperture  and  the  signal  ring  Rj  and  wedge  W.  values  are  generally 
of  the  same  order  of  magnitude.  As  the  aperture  is  made  smaller,  we  see  that  the  S/N 
decreases,  the  positive  slope  of  the  signal  detector  elements  Rj  and  Wj  increases,  and 
the  negative  slope  of  the  adjacent  detector  elements,  R:j,  RJ+i,  and  WJ+).  decreases. 
For  the  smallest  sampling  aperture  at  this  focal  length.  0.1  mm,  and  S/N  values  for  the 
wedges  are  positive  throughout:  whereas,  the  S/N  values  for  the  rings  are  less  than  one 
for  the  lower  spatial  frequencies. 

For  a  transform  lens  focal  length  of  100.0  mm,  the  same  generalizations  are 
valid  with  the  additional  statem-mt  that  the  S/N  values  are  higher.  This  is  readily  seen 
by  comparing  figures  24  and  21 . 


At  this  time,  we  will  consider  the  details  of  curve  configuration  and  curve 
differences,  Graphical  and  tabular  summaries  will  be  introduced  to  aid  in  this  analysis. 


Fhe  similarity  of  figures  2 ;  and  25  is  important.  We  note  that  for  figure  21 . 
when  f  =  177.8  mm  and  1)  =  0.5  mm,  the  relative  aperture,  given  by  f*  =  f  d.  is 
f/555.6:  and  for  figure  24.  when  I  =  100.0  nun  ami  I)  =  0.25  mm.  the  relative 
aperture  is  17400.0.  As  shown  in  table  1.  the  relative  aperture  provides  the  scaling  func¬ 
tion  for  the  relative  si/e  of  the  Bessel  function  in  the  spatial  frequency  domain.  The 
Bessel  function  zero  crossings  pertaining  to  Data  Set  1  can  be  abstracted  from  table  1  and 
shown  in  table  6  as  diameters  between  Bessel  function  zeros.  Considering  that  85.6 
percent  of  the  energy  of  the  Bessel  function  is  contained  in  the  central  maximum,  i.e.  the 
diameter  between  the  first  zero  crossing,  that  only  7.1  percent  of  the  energy  is  between 
the  first  and  second  zero  crossings,  and  that  2.8  percent  and  1.5  percent  of  the  total 
energy  occurs  between  the  second  to  third  and  the  third  to  fourth  zero  crossings,  res¬ 
pectively.  one  can  use  the  detector  ring  and  wedge  Fmensional  data  in  table  2.  with 
the  Bessel  t unction  zero  crossing  dimensions  to  assist  m  understanding  the  ring  and  wedge 
data  in  figures  21  through  27. 
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t  =  1 77.8  mm 

D  =  0.250  mm 

WRD-6420  DETECTOR 


RING  NUMBER 

'  W  ’  '  T  1  r  ▼"*  *  ■  *  •  »  ■  I  T1  »'  '  T*  r-  I  — 

7  8  *  10  >2  14  I*  18  20  22  24  2*  28  30  32 

FKJlIRh  22,  ROSA  Simulation  Results.  B 
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f  =  177.8  mm 

U  -  0.1  mm 

WRD— 6420  DETECTOR 


RING  NUMBER 


7  8  9  10  12  14  16  18  20  22  24  26  28  30  32 

FIGURE  23.  ROSA  Simulation  Results.  C 
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10  2 


|  f  =  100.0  mm 

10  31  1>  =  0.075  mm 

WRD-6420  DETECTOR 


7  8  9  10  12  14  16  18  20  22  24  26  28  50  52 


FIGURE  27.  ROSA  Simulation  Results.  G 
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Anotlier  useful  type  of  data  for  interpretation  of  Data  Set  1  is  the  location 
of  the  Ronchi  Grating  deltas  for  the  mid-ring  spatial  frequencies  associated  with  the 
WRD-6420  detector  when  the  transform  lens  focal  lengths  are  specified.  Table  7  re¬ 
presents  data  of  this  type  for  a  177.8  mm  transform  lens  focal  length,  and  table  8  has 
similar  data  for  a  100.0  mm  lens. 


The  next  data  analysis  aids  to  be  introduced  are  the  generalized  ring  and 
wedge  models  shown  in  figure  28.  The  ring  model  shows  the  Bessel  functions,  represented 
in  plan  view  by  the  four  zero  crossing  circles,  located  at  the  mid-ring  radius,  and  each 
successive  ring  has  a  greater  radius  difference,  as  shown  in  table  2  and  figures  17  and  18. 
The  scale  of  the  Bessel  function  is  arbitary  and  can  be  adjusted  to  satisfy  our  interpre¬ 
tation  of  the  mod  nd  to  fit  the  curve  being  interpreted.  The  wedge  model  of  figure  28 
represents  the  Bessel  functions  positioned  for  a  Ronchi  grating  spatial  frequency  of 
30.34  line-pairs  per  millimeter  when  the  177.8  mm  transform  lens  is  employed.  This 
corresponds  to  the  data  for  Ring  13  in  table  7.  This  wedge  model  has  3  general  cases 
and  18  special  cases,  shown  in  figures  29  and  30.  respectively.  The  general  cases  apply 
where  all  of  the  Bessel  function  energy  within  the  fourth  zero  radius  is  incident  on  the 
detector  wedges,  f  or  the  first  general  case,  the  Bessel  function  is  scaled  so  that  the  total 
significant  energy  is  detected  by  W;  and  the  energy  in  Wj.j  and  WJ+|  is  considered 
insignificant.  In  general  case  two,  all  of  the  energy  is  contained  in  wedges  Wj.,,  Wj. 
and  Wj+1  ;  and  for  general  case  three,  the  energy  is  incident  on  other  adjacent  wedges. 
The  special  cases  apply  when  portions  of  the  scaled  Bessel  function  exist  outside  of 
the  detector  wedges.  In  figure  30,  the  special  cases  of  the  bottom  row  apply  when  the 
lower  radius  of  the  detector  wedges  (WLR)  is  greater  than  the  axial  distance  of  the 
delta  harmonic  r(,  less  the  fourth  zero  radius  of  the  Bessel  function  r,.  Likewise,  the 
special  cases  in  the  top  and  middle  rows  of  figure  30  apply  when  portions  of  the  fourth 
zero  circle  are  beyond  the  higher  radius  of  the  detector  wedges  (WHR).  Table  9  represents 
a  summary  of  the  occurrences  of  special  cases  for  selected  simulation  with  the  WRD-6420 
detector,  which  will  be  useful  in  analyzing  curves  for  detector  signals. 


f  igure  21  can  be  analyzed  as  a  simple  example  of  a  detailed  curve  analysis. 
The  remaining  curves  in  Data  Sets  1 A  and  IB  can  be  analyzed  in  like  fashion.  First,  note 
that  the  scaled  Bessel  function  at  the  origin  has  a  fourth  zero  crossing  radius  less  than  the 
least  radius  of  Ring  7  (tables  6  and  2).  Next,  note  that  with  the  first  harmonic  positioned 
at  the  mid-radius  of  Ring  7,  there  will  be  no  signal  contribution  on  Ring  7  owing  to  the 
Bessel  function  positioned  at  the  third-order  harmonic  (tables  2.  (>  and  7).  This  lack  of 
signal  is  easily  seen  when  the  third-order  Bessel  function  is  positioned  at  an  axial  dis¬ 
tance  of  3.855  mm  and  when  the  radius  of  fourth  zero  crossing  for  f  =  177.8  mm 
and  I)  =  0.5  mm  is  0.954  mm.  When  subtracted,  the  resulting  radial  distance 
<2.901  mm)  is  greater  than  the  maximum  ring  radius  (1.395  mm)  for  Ring  7.  Thus. 
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TAB_E  7.  Summary  of  Ronclii  Grating  Deltas  for  177.8  mm  Transform 


u 

4-> 

•r— 

•r— 

C 

CO 

o 

na 

c 

E 

4— 5 

<D 

i_ 

r— 

4-> 

0) 

C 

:r 

Q 

► — i 

o 

r  i 


i 

X. 

£ 


O 

\> 


S'. 


co 

o  o 

f-H  O  CO  CO 

rH  O  T— t  CO 

O  <3-  r-. 

'  X  rH  rH 

5- 

O 


4-> 

rH 

^3"  C"- 

CM  CO 

U 

as  o 

to  rH 

CO  C"- 

- - X 

aj 

o 

to  CM 

O  t-H 

E 

4-> 

•  * 

• 

E 

a) 

o 

t-H 

O 

- ’ 

a 

rH  rH 

*— H  CM 

CO 

-o 

u 

c 

CNJ 

•X  00 

r^. 

as  cm 

•r- 

o 

r-x  o 

os  tn 

oo  as 

t-H 

c 

>) 

o 

as  o 

CM  to 

oo  o 

o 

QJ 

•  • 

> 

E 

CQ 

o 

CO  r-H 

m  m 

CO  r-H 

i- 

' — t 

r-H  r-H 

rH  CM 

«3 

co 

zn 

o 

•r- 

*3-  CO 

O  CM 

to  CM 

to  o 

o  o 

♦r~ 

C 

LO  O 

cm  as 

as  i — i 

CO  oo 

as  co 

rx 

-C 

o 

o 

*3-  co 

CM 

O  O 

O  CM 

LO  os 

U 

E 

• 

*  • 

•  • 

•  • 

•  • 

c 

i- 

o 

to  r-x. 

OS  i — t 

CO  LO 

h-  as 

rH  CO 

o 

OJ 

, — < 

rH  rH  CM  CM 

C£ 

3: 

<4- 

03 

» — < 

to  os 

<3-  r-x 

r-H  OS 

CO  Cx. 

rH  CM 

» — <  to 

r^- 

r*x 

o 

00 

o 

t— H 

t — i 

r- 

as  »— h 

» — < 

rH 

<3- 

to 

O 

i- 

CO  rH 

LO  co 

as  cm 

LO  rH 

tn  to  ^  co 

as  to 

o 

rH 

as 

co 

S3* 

i — 1 

LO 

LO 

CM  to 

<3- 

Cx. 

o 

00 

U3N  030 

cm  to 

as  ro 

co 

t-H  00  LO 

oo 

co 

CO 

to  CO 

— ♦ 

to 

as 

o 

CO 

o 

oo  <3* 

LO  tO 

r-  as 

«H 

CM  ^ 

LO  r~x 

os 

O 

cm  ^3- 

o 

CO 

o 

CM 

to 

as 

« — ( 

CO 

to 

c 

a> 

i — 1  i — 1 

- - 1  tH 

< — »  t-H 

rH 

CM 

CM 

CM 

CM 

CM 

oo 

OO 

CO  CO 

CO 

*3- 

^3- 

X 

u 

4-> 

c 

CO 

03 

rH 

to  o 

co  cm  to 

as  to 

*7T  ^ 

cm 

as 

as 

o  to 

O  rx. 

r-x 

as 

oo  r- 

r-x 

CO 

CM 

•r- 

4~> 

rH  O 

co  00 

as 

rH  O 

r-H  LO 

i — l  OS 

OS  CM 

to 

OO 

oo 

co  oo 

rH 

*H  OO 

LO 

LD 

as 

o 

CO 

T— « 

CM  UO 

CO  CM 

to  O 

*3-  CO 

co  rx 

CM  CO 

CO 

as 

LO 

t-H 

r-x. 

*73- 

rH 

CO 

LO  CM 

o 

co 

to 

r- 

Q 

o 

rH  r-H 

t-H  CM 

CM  00 

CO  co 

*73-  *73- 

LO  LO 

to 

to 

CO 

co 

as 

o 

o 

r — i  CM 

OO 

co 

lO 

03 

1 — t 

rH 

rH  r — i 

f— 1 

•t— 

"O 

X 

a> 

c 

i — 

o 

o 

O  to 

o  o 

0 

0 

0 

0 

0 

0 

o  o 

0 

0 

o 

o 

o 

o 

O 

o 

o 

o 

0 

0 

o 

o 

O 

o 

s- 

CM 

C-J  o 


as 

c 

f _ 

>s 

U 

CM  ^ 

Ox-  co 

tO  CM 

as  r — 

rH 

co  tn 

rH 

r-x 

CO 

o 

oo 

CM 

to 

as 

< _ ( 

to 

Cx. 

CO 

Os 

•r— 

03 

C 

^  O 

co  as 

cm  r-x. 

CO 

CM 

CO  to 

O 

to 

to 

as 

rj- 

o  cc 

as 

rH 

tn 

CM 

o 

o 

CM 

to 

Cd 

•  r— 

03 

•  • 

•  • 

•  • 

• 

• 

•  » 

•  • 

• 

• 

• 

« 

• 

• 

1 

4-J 

3 

rH  ^ 

to  OS 

co  to 

o 

00  CM 

N  rH 

to 

rH 

to 

CM 

CO  co 

os 

to 

CM 

Os 

to 

oo 

o 

r-x 

*o 

<T3 

an 

rH  rH 

rH  rH 

CM  CM 

CO 

CO 

CO  ^ 

tn 

tn 

to 

to 

f^x 

r^x  co 

oo 

as 

O 

O 

rH 

CM 

co 

ro 

£ 

Q. 

to 

03 

i- 

Li_ 

rH 

rH 

r  1 

rH 

r-H 

rH 

as  « 
c  o 

r-x  CO  OS  o  rH  CM 

oo  in  to 

fXx  00 

OS 

o 

rH  CM 

CO 

*73- 

tn 

to 

f^X 

00 

os 

O 

t-H  CM 

•f-  2T 

Cd 

rH  rH  rH 

rH  rH  rH  rH 

rH  rH 

t-H 

CM 

CM  CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

os 

CO  CO 

65 


TABLE  8.  Summary  of  Ronchi  (.rating  Deltas  l  or  1 00.0  nun  Transform 
Lens  Focal  Length  (Detector  WRD-6420). 
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LET:  r  *'■  axial  distance  of  i,h  deha  harmonic. 

P 

r.  =  radius  of  inscribed  circle  in  wedge  W  centered  at  r. 

I  -  L  J  P 

r2  =  radius  of  Bessel  function  4'  rero. 

r3  =  radius  of  incribed  circle  in  wedges  W.  t  +  Wj  +  Wj+1.  centered  at  rp. 


GENERAL  CASE  #1  -  When  r2  ^  r,  . 

Energy  Ej  on  wedge  Wj  due  to  Bessel  function  centered  at  rp  is 
Ej  *  1,  |1  -  J,2,  ftr  rjD/X  f)|  =  \,2 

where  I,  is  the  intensity  of  the  ilh  delta  harmonic. 

Energy  Ejh  on  Wedge  WJlf,  is  considered  insignificant. 


GENERAL  CASE  #2-  When  r,  <  r2  rr 
Energy  E(  on  wedge  Wj  is 

E.  *  '  '.p' 

where  P’  is  percentage  of  Bessel  function  energy  in  sectors  outside  of  Wj. 
Energy  E)h  on  wedge  WJ(1 

t,h  55  ‘i  V'l2 


GENERAL  CASE  #3  -  When  r2  >  r3 
As  in  General  Case  #2  Energy  on  W(  is 

h  *  '.2  -  «iPp’ 

Energy  Ej,,  on  wedge  Wj4|  is 

E|* ,  *  1,  <P'-P")/2 

w  here  P "  is  percentage  of  Bessel  function  energy  in  sectors  outside  of 
(W(,,  +  W,  ♦  W,,,). 


lIGlIRli  29.  Three  Ueneral  Cases  for  Estimating  Bessel  I  unction  Energy 
Incident  on  a  ROSA  Wedge  Detector  Element. 
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FIGURE  30.  Eighteen  Special  Cases  lor  Estimating  Bessel  Function  Energy 
Incident  on  a  ROSA  Wedge  Detector  Element. 


the  ring  detector  signals  for  conditions  of  figure  21  are  due  to  a  single  Bessel  function 
positioned  at  the  mid-ring  radius.  The  IT  curve  is  directly  relatable  to  the  plots  of  ring 
widths  in  figure  17  and  18,  and  the  shapes  of  R-_,  and  R|+ ,  curves  can  also  be  related 
to  these  plots. 


The  wedge  signals  for  figure  21  can  be  analyzed  from  table  7,  which 
defines  the  axial  distances  and  intensities  of  Ronehi  deltas,  from  figures  29  and  30, 
which  describe  the  general  and  special  cases  for  estimating  wedge  signals,  and  from  table 
9,  which  summarizes  the  occurrences  of  the  special  cases.  It  will  be  noted  that  all  W( 
curves  in  figures  21  through  27  have  a  discontinuity  between  Ring  lb  and  Ring  17.  This 
discontinuity  is  due  to  the  transition  from  the  occurrence  of  the  first  and  third  harmonics 
on  the  detector  for  the  Ring  16  data  to  the  first  harmonic  on  the  detector  when  the 
third  harmonic  is  off  of  the  detector  for  the  Ring  17  data,  as  seen  from  tables  7  and  9. 
The  slope  changes  in  the  Wj+(  curves  can  be  related  to  the  diameter  of  the  Bessel 
function  zeros  in  table  6,  the  wedge  width  data  in  table  2,  and  the  slope  changes  in 
figure  10. 


Table  10  summarizes  the  number  of  detector  elements  with  S/N  greater 
than  1.0  for  the  data  in  Data  Set  1A.  It  should  be  noted  that  when  the  f#  is  less  than 
about  1 .000,  all  26  rings  and  32  wedges  have  S/N  greater  than  1 .0,  and  the  S/N  improves 
inversely  with  f#. 


2.  Data  Set  IB -Detector  WRD  6400  Simulations.  Table  5  lists  the 
simulation  conditions  under  which  Data  Set  IB  was  obtained.  The  essential  differences 
between  these  data  and  those  of  Data  Set  1A  were  the  detectors,  which  have  been  des¬ 
cribed  in  the  section  Geometry  of  ROSA  Detectors.  Data  Set  IB  simulation  results 
obtained  with  the  WRD-6400  detector  are  presented  in  figures  31  through  36.  and  a 
quick  comparison  with  their  counterparts  in  Data  Set  1A  will  show  significant  differ¬ 
ences.  The  simulation  results  of  Data  Set  IB  are  generally  less  well-behaved  and  have 
lower  S/N  ratios.  Details  of  the  plots  can  be  explained  by  using  the  plots  of  mid-ring 
radii  for  the  two  detectors  in  figures  15  and  16,  the  plots  of  the  ring  widths  in  figures 
I  7  and  18.  the  dimensions  of  the  two  detectors  in  tables  2  and  3.  and  the  summaries  of 
Ronehi  grating  deltas  for  the  177.8  and  100.0  mm  transform  lens  focal  lengths  shown 
in  tables  1  1  and  1  2.  respectively . 
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TABLE  10.  Detector  Elements  with  S/N  >  1  In  Data  Set  1  A. 


f 

D 

Fig.  21 

177.8  mm 

0.5  mm 

22 

177.8 

0.25 

2: 

177.8 

0.10 

24 

100.0 

0.5 

25 

100.0 

0.25 

26 

100.0 

0.10 

27 

100.0 

0.075 

Number  of  Detector 
Elements  with 


S/N>1 

f/D 

Ri  ngs 

Wedges 

355 

26 

32 

711 

26 

32 

1778 

21 

32 

200 

26 

32 

400 

26 

32 

1000 

25 

32 

1333 

23 

32 
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I  IGURb  31.  ROSA  Simulation  Results.  II 
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TABLL  1 1 .  Summary  of  Ronchi  Grating  Deltas  For  1  77.8  mm  Transform 
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T  \BLF.  12.  Summary  of  Ronchi  Grating  Deltas  E'or  100.0  mm  Transform 

Len.  T  .al  Length  (Detector  WRD-6400). 


The  simulation  data  for  the  WRD-6400  detector  becomes  particularly  ill- 
behaved  for  both  transform  lens  focal  lengths  when  the  sampling  aperture  has  a  diameter 
of  0.1  mm  or  less.  An  example  of  the  manner  in  which  a  S/N  value  less  than  1.0  mm 
is  obtained  for  the  f  =  100.0  mm,  1)  =  0.1  min  data  shown  in  figure  35  is  illustrated 
in  figure  37  where  the  sections  through  the  zero  and  first-order  harmonics  have  been 
plotted  for  the  first  harmonic  centered  on  Ring  12.  If  we  algebraically  add  the  two 
curves,  square  the  results,  and  then  integrate  over  the  Ring  12  and  Ring  13  areas,  we 
find  that  the  R  /R  ,  is  less  than  one. 


It  will  also  be  noted  that  the  wedge  data  have  lower  S/N  ratios  for  the 
WRD-6400  detector  because  of  the  smaller  wedge  angles  when  compared  to  those  of  the 
WRD-6420  detector. 


Table  13  summarizes  the  number  of  detector  elements  with  S/N  greater 
than  1.0  for  Data  Set  IB.  When  this  data  is  compared  with  the  similar  data  in  table  10 
for  the  WRD-6420  detector,  we  see  more  rings  with  S/N  greater  than  1.0,  but  it  must 
be  remembered  that  the  S/N  are  in  general  lower,  especially  for  the  lower  and  middle 
detector  rings,  as  shown  in  figures  38  and  30  for  f  =  100.0  mm  and  D  =  0.1  mm. 


3.  Data  Set  2A  -  Ronchi  Grating  Deformation.  Table  5  lists  the  simula¬ 
tion  conditions  under  which  Data  Set  2A  was  obtained.  The  purpose  of  this  data  set  was 
to  examine  the  change  in  S/N  when  the  Ronchi  grating  was  subjected  to  simple  linear 
deformations.  These  linear  deformations  were  obtained  by  multiplying  the  mid-ring 
spatial  frequency  by  a  factor  that  had  a  value  between  0.06  and  1.04.  Note  that  this 
is  not  the  equivalent  to  selecting  a  deformation  and  then  computing  the  S/N  values  for 
each  ring  detector  element,  l  or  example,  a  4  percent  deformation  having  its  first  har¬ 
monic  centered  on  Ring  7  for  f  =  100.0  mm  with  the  WRD-6420  detector  is  only 
0.051  mm;  whereas,  the  same  percent  of  deformation  for  the  Ring  32  spatial  frequency 
would  amount  to  a  deformation  of  0.620  mm.  In  either  case,  however,  the  effect  of 
the  deformation  is  to  move  the  position  of  the  first  harmonic  grating  function  from  its 
mid-ring  position.  Figures  38  and  39  plot  the  mid-ring  S/N  values  for  Rj/Rj.,  and 
R,/RJ+ 1  for  the  WRD-6420  and  the  WRD-6400  detector  respectively  when  f  =  1 00.0  mm 
and  D  =  0,10  mm.  Table  14  lists  the  S/N  values  for  grating  deformations  associated 
with  the  WRD-6420  detector,  figure  38  plots  S/N  values  for  no  deformation  (f  =  1.0) 
while  figure  40  plots  S/N  for  selected  rings  from  table  14.  Table  15  is  similar  to  that 
of  table  14.  except  that  D  =  0.075  mm  instead  of  0.10  mm. 
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amplitude 


TABLE  13.  Detector  Elements  with  S/N  > 


f 


Fig.  31 

177.8  mm 

32 

177.8  mm 

33 

177.8  mm 

34 

100.0 

35 

100.0 

36 

100.0 

D 

f/D 

0.5  mm 

355 

0.25  mm 

711 

0.10  mm 

1778 

0.25  mm 

400 

0.10  mm 

1000 

0.075  mm 

1333 

1  In  Data  Set  IB. 


Number  of 

Detector 

Elements 

with 

S/N>  1 

Rings 

Wedc 

27 

32 

27 

32 

26 

32 

27 

32 

26 

32 

26 

32 

i 

1 
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Note:  S/N  lor  K/K,.,  anil  R,/R)f|  When  f  =  100.0  nun. 
anil  WRD-M20  1)1  I  I  C 'TOR. 


I'KiURK  38.  Signal-to-Noise  I’lolv  A 


1 00  I 
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6  7  8  9  10  12 


14  16 


18  20  22  24  26  28  30  32 


Note:  S/N  for  / RJ . ,  and  R;/RJ+J  When  f  =  100.0  mm,  1)  =  0.1  mm, 
and  WRD-6400  DETECTOR. 


FIGURE  39.  Signai-to-Noise  Plots.  B 


85 


I  —*000000  00000 


1  M  cvj  m  (\j  n  n  n  pi  n? 


^  «y  in  in  in  i/i  in 


1  cvj  c\j  n  m  uiio  r 


.o>OMM«ji/)inioioioio  m  10  o  ud  m 


.  m  n  m  tt  «d-  *3-  oeoroo  c\j  <\j  cm  c\j  cm 


•— «  (\j  m  i\j  fvj  m-  in  r 


f  in  ns  co  O'  o  rHW  nmnn  r--~  co  0  o  «  cm 
cm  cm  cnj  cm  cm  cm  cm  (\j  cm  nj  <n  n  ,0 


.  m  id  m  ui  o  <?  *3-  <3-  < 


riniDrscoo^cccOWcocco^cOCCCO 


uji/uONa'QiMnuMOr 


cm  r-  cm  f..i  m  m  >:f  «.t  in  10  10  ncO  d  . 


-  CO  <T>  O  O  ■ — ■  < — •  cm  r 


<-  «ij  «:j-  tf  inin  un  in  in  in  1 


IPJCMIMCMCMIMCMIMIMCMIMCM  CM  CM  f 


.-..-<,-.0  0  0 


•  10  in  if>  in  »t  co  r vj  0.1  . 


.  —*  —«  o  o  o  o  o  oc 


0.96  0.97  0.98  0.99  1.00  1.01  1.02 


.03  1 .04 


FACTOR 

FIGURE  40.  Plots  of  Selected  Data  From  Tabic  14. 

87 


3 

O  CO 

cm  r- 

VO  O  00 

3 

•31 

CM  O 

r-. 

3  cn 

VO 

r-H 

LO  O 

*3- 

cn  co 

CO  co 

r^. 

ro 

o 

o 

r-«. 

CO  co 

ro  t-H  O 

i  i 

o 

o  cn 

cn 

cn 

CO  00 

C-. 

vo  vo 

to 

LO 

3 

3 

- 

fl  1 

O 

O  r-H 

—•  •“*  f_* 

-1 

— 

t—i 

t“H  O 

0 

0 

O  O 

O 

00 

0 

O 

O 

0 

<\j 

LO 

o 

CTl  3 

co  ro 

CO 

ro 

VO 

r-H  CM 

r»- 

00 

cn  r-H 

r-H 

CM  O 

CO  VO 

CO 

t-H 

o 

O  CO 

CO  u*> 

lo  co  ro 

LO 

LO  VO 

VO 

r>- 

r-» 

r- 

r>»  00 

co 

coco 

1 — . 

p— 

r-» 

r~< 

f— i 

o 

O  > — 1 

- 

r-H  r-H 

— 

r-H  r-H 

— 

— 

* — i  T—i 

r— H  r-H 

r-H 

■ — i 

— 

t-H 

*  ■  4 

CO 

CTl 

ro  ro 

rH  VO  CTl 

CO  VO  CO 

CO 

O  CO  CM 

CM 

cm  3 

r-»  lo 

lO 

«3" 

CM 

o 

O  CO 

CT  VO 

3  3 

vo 

r>. 

co  cn 

o 

CM  CM 

<3- 

LO 

vo  r-. 

co 

<n  0 

CM 

ro  *3- 

- 

rH 

o 

O  t-H 

— 

-■ 

-H 

CM 

CM  CM 

CM 

CM 

CM  CM 

CM 

CM  co 

ro 

CO 

CO 

ro 

o 

ro 

r-»  ro 

■d-  co 

co 

CM  CM 

r-H 

LO  LO 

CO 

CO  VO  CO 

00 

to  co 

3 

CO 

cn 

CM 

o 

*-* 

CTl 

CT>  r-*. 

co  cn  r-H 

cn 

o 

cm  3 

vO 

CO  O 

CO 

LO 

00  t-H 

CO  »— 1 

LO 

cn 

CM 

r--. 

O 

O  i 

r-  H  r  4  r-H 

r— H 

CM 

CM  CM 

CM 

CM  CO 

CO 

CO 

ro  *3- 

*3- 

3  LO 

LO 

LO 

VO 

vo 

cr» 

ro 

O  3 

CTi  VO  VO 

co  co 

VO  CO 

O 

o  cn 

1— H 

CM 

r-.  co 

cn 

00  0 

0 

VO 

co 

O 

cr> 

rH 

<n 

O  CO 

cti  r^.  co 

—* 

3 

vo  cn 

CO 

r-  o 

vo 

*— • 

vo  CO 

cn 

r-^  lo 

ro 

cn 

cn 

o 

o 

_  __ 

r-H  r-H  rH 

CM 

CM 

CM  CM 

ro 

ro  *3- 

3 

LO 

LO  VO 

vo 

r-  co 

cn 

O 

0 

t-H 

t-H 

r-H 

CO 

in 

C\J 

3  in 

UO  3"  CTl 

O  to 

O  CO 

cn 

3  CO 

vo 

CM 

LO  CM 

LO 

0  vo 

r— H 

CM 

to  00 

CTl 

r-H 

o 

O  cn 

rH  CTl  O 

vo  CO  csj  vo 

•  < 

co  3 

ro 

CM 

r—i  CM 

CM 

«3"  CO 

CO 

CM 

cn 

VO 

o 

—t 

CM  r-H  CM 

cm  oj  ro  ro 

3  uo 

vo 

CO  CTl 

0 

r-H  CM 

CO 

vt 

3 

LO 

t-H 

r-Hr-H 

r-H 

t-H 

■“* 

VO 

CM  * — i 

cun  vr 

r-H 

VO 

r-  vo 

CM 

*3-  vo  CO 

10 

*3-  O 

r-H 

<0  r- 

10 

CM 

CO 

CM 

CTl 

CM 

r-H  CM 

lo  ro  cO 

co 

cn  vo  'O  r^. 

O  CM 

vo 

3  CO  CO 

VO  r-H 

r-. 

CO 

CO 

LO 

o 

r-H 

r-H 

t” H  CM 

CM  CM  CM 

co 

co 

LO 

vo  co  cn 

O 

r—i 

CM  CM 

CM 

CM  CM 

T—i 

r-H 

0 

O 

t — i 

r— H 

>—i  t—i 

*“H 

t — t  t — i 

r—i 

r-H 

t-H 

CTl  ♦ 

C  O 

co 

CT>  o 

r-H  CM  ro 

3 

uo 

vo  r-* 

co 

cn  0 

1 — 1 

CM 

CO  3 

LO 

vo 

00 

cn 

O  t-H 

•r-  Z 

t-H 

t—H  r-H  r-H 

r-H 

*— H 

f—t  r-H 

•— * 

«— »  CM 

CM 

CM 

Ol  CM  CM 

CM  CM 

CM 

CM 

ro 

CO  . 

0  0 

r- 

CM  r-H 

vo 

CM  CM 

CO 

O 

CM  CO 

CM 

cn 

r-  CO 

< — 1 

CM  r-H 

VO  CO 

m 

0  VO 

0 

CM  O 

cn  co 

CO  <03 

cn 

t-H  CM 

ro 

U1 

cn  3 

CM 

-3- 

ro  co 

00 

3  r-*. 

O  3 

co  3  cn 

• — *  CM 

CM  CM 

r-H 

cm  ro 

CO 

LO  LO 

vo 

c- 

CO  O  CM 

CO  *3  *3" 

3 

3  ro 

ro 

CM 

•-H 

t-H  O 

r-H 

r—i  r—i 

—• 

t-H  r—i 

r-H 

t—H 

f—H  r—i 

CM 

CM  VO 

CO  vo 

LO  *3-  3 

CM  CO 

t-H 

ro 

vo  CO 

VO 

O 

O  r-H 

CO  3 

t-H 

00 

CM 

cm  r-*. 

O 

r-H  cn 

3  VO 

VO  CTl  N 

O  CM  vo 

vo 

cm  cn 

CTl 

CTi 

CTl  r-H 

3 

vo  r^. 

CTl 

00 

VO 

co  r- 

r-H 

r-H  t-H 

CM  CM 

r-H 

r—i  CM 

CO 

CO  CO 

3 

*3 

LO  LO 

VO 

r». 

00  0 

t—H 

cm  ro 

3 

in 

vo 

r->  r-. 

*“• 

r—i 

t—H  M 

t—H 

t-H 

r—i  r — i 

r-H 

cn  <3- 

3  CM 

LO 

r-»  vo 

co  vo  cn 

*3 

VO 

CO  *3 

t-H 

3 

O  O 

O 

cm  cn 

LT> 

0 

co  in  cm 

O 

0  CO 

CO  <3- 

U1 

r*»  3 

vo 

f--  CTl 

CO 

VO 

0  -3- 

O 

LO 

t-H  CO 

VO 

3  t-H 

t—H 

t-H 

CTl 

CTi  CO 

T-i 

t-H  r—i 

CM  CM 

r-H  CM 

CM 

CM  CM 

CO 

CO 

3  3 

LO 

LO 

vo  vo 

CO  CTi 

O 

t-H 

t—H 

CM  CO 

t— 1 

r-H 

1 — 1 

t-H  r-H 

O 

LO  CO 

*3  CO 

LO 

ro  vo 

CM  r-.  r-H 

CM 

0 

O  CM 

CM 

r-- 

CM  «3 

t—H 

r-.  0 

CM 

CO 

CM  3  CM 

O 

O  P's. 

CM  CM 

*3- 

VO  LO 

CO 

CO  LO 

r-. 

<n 

r-H  ro 

VO 

co 

t-H  3- 

00 

t-H  m 

cn 

ro 

t-H  m 

r-H 

CM  CM 

CM 

CM  CM 

CM 

CM 

ro  ro 

CO 

ro 

3-  3- 

3 

vn  in 

in 

vo 

vo 

cn 

CM  CM 

LO  VO 

r-. 

CTl  CTl 

CO  LO  CM  3 

r-H 

t-H  0 

CO 

CM 

O  t—H 

3 

3  CM 

co 

co 

0 

cm 

cn 

0  r- 

t-H  r-H 

CO  <3-  CTl 

O 

O  fH 

CM 

ro 

3  LO 

vo 

r-» 

co  cn 

O 

r-H  CM 

ro 

3 

in 

in  vo 

0 

r-H  t-H 

CM  CM 

*— H 

^H  t—H 

CM 

OJ  CM 

CM 

CM 

CM  CM 

CM 

CM 

CM  CM 

ro 

ro  ro 

ro 

CO 

co 

ro  co 

00 

CTi  VO 

VO  3 

co  n  cn 

3  cn  3 

vo  co  cn 

ro 

3 

ro  o"' 

3 

ro  O 

CTl 

vo 

CM  00  co 

cn 

cn  vo 

O  O 

CM 

ro  c-. 

co  r-  N  CO  CO 

CO  00  CTl 

cn 

cn  cn 

<n 

CTl  CTl 

co  00  CO  n 

0 

O  »-H 

CM  CM 

t-H 

t-H  r-H 

- 

r—i  t—i 

t-H 

r-H 

t — H  t—i 

t—H 

t-H 

t—H  r—i 

r-H 

r-H  r-H 

r-H 

f— H 

- 

t—H  r-H 

vo 

3  3 

O  3 

ro  vo  to 

10 

n  cn  n 

^H 

vo  0 

VO 

O 

ro  r** 

t-H 

3  00 

CM 

vo 

O  LO  O 

cn 

CTl  VO 

cn  co 

•-* 

t-H  3 

CO  CM  CM 

CM 

t-H  T—H 

0 

O  CTi  00 

co  r-%  vo 

VO 

m 

VO  3  3 

0 

O  r-H 

r—i  r-H 

r-H  r-H 

t-H  t-H 

t-H 

r-H 

t-H  i-H 

t-H 

^H 

O  O 

0 

OO 

0 

O 

O 

O  O 

0 

r-.  co 

cn  0 

t-H  CM  fO 

lO  VO  N* 

CO 

cn  0 

CM 

ro  3 

10 

VO 

00 

cn 

O 

t-H  CM 

2 

t-H 

t—H 

r-H  r-H 

r-H 

r— H  t-H 

t-H 

t—H 

r-H  CM 

CM 

CM 

CM  CM 

CM 

CM  CM 

CM 

CM 

CO 

CO  CO 

88 


4.  Data  Set  2B  -  Ron  chi  Grating  Rotation.  This  data  set,  as  shown  in 
table  16,  was  obtained  to  investigate  S/N  values  when  the  Ronchi  gratings  were  rotated 
from  the  wedge  bisector.  The  f  and  D  conditions  under  which  the  data  were  obtained 
are  the  same  as  figure  38  and  table  14.  The  results,  as  expected,  show  a  decrease  to  a 
S/N  value  of  1.0  when  the  grating  is  rotated  2.3°.  neglecting  a  small  but  finite  seperator 
between  wedges. 


5.  Data  Set  3A  -  Hqual  Ring  Width  Detector.  In  section  4,  the  simulation 
data  indicated  response  differences  for  the  two  detectors.  This  motivated  an  attempt 
to  redesign  the  ROSA  detector  to  perform  better.  Data  Set  3 A  represents  simulation 
data  horn  a  detector  with  equal  ring  spacing  when  D  =  0.1  mm  and  the  transform  lens 
focal  length  were  100.0,  150.0,  and  200.0  mm. 


A  summary  for  a  ROSA  type,  equal  ring  width  detector  is  tabulated  in 
table  17.  This  detector  had  32  rings  and  32  wedges,  with  the  mid-ring  radii  0.5  mm 
apart  and  with  0.0254  mm  (1.0  mil)  allowance  for  detector  element  sperators.  Table 
17  gives  the  mid-ring  spatial  frequencies  for  f  =  100.0  mm.  f  =  150.0  mm.  and 
f  =  200.0  mm  and  the  Ronchi  grating  delta  positions  when  the  first  harm*  lie  was 
positioned  at  the  mid-ring.  These  detector  data  can  be  compared  with  similar  data  in 
tables  2.  3,  7.  8,  11.  and  12. 


Simulation  data  for  the  equal  ring  width  detector  are  plotted  in  figures  41 
through  43  for  D  =  0.1  mm  and  f  =  100.0  mm,  f  =  150.0  mm,  and  f  =  200.0  mm, 
respectively.  For  a  quick  appraisal  of  the  results,  compare  figure  41  with  figure  26  and 
figure  35,  and  note  that  although  the  absolute  response  is  reduced  for  the  rings  and 
wedges,  the  S/N  is  increased  and  almost  constant  for  the  rings.  Additionally,  note  that 
rings  3  through  32  and  wedges  3  through  32  all  have  S/N  greater  than  1.0.  Figures  42 
and  43  indicate  that  S/N  decreases  with  increased  transform  lens  focal  length,  and  the 
noise  seen  in  figure  40,  for  the  lower  rings,  becomes  more  sever  as  I  increases.  Table  IS 
summarizes  the  number  of  detector  elements  with  S/N  greater  than  1.0  for  Data  Set  3A. 


6.  Data  Set  3 H  -  Modified  Fqu.i!  Ring  Width  Detector.  Although  there 
were  noisy  variations  in  the  simulation  data  for  Data  Set  3A  for  the  lower  rings,  a  brief 
attempt  was  made  to  adjust  ring  width  to  minimize  this  noise.  However,  these  attempts 
were  unsatisfactory,  and  the  data  appears  in  appendix  B  as  indicated  in  table  5. 


TABLE  16.  S/N  For  Grating  Rotations  From  Wedge  Bisector,  f  =  100.0  mm. 

D  =  0.1  mm.  (Detector  WRD-6420). 
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V.  PRELIMINARY  OPTICAL  SUBSYSTEM  DESIGN 
FOR  GRATING  MAP  SYSTEMS 


In  the  preceding  section,  we  determined  that  Ronchi  gratings  having  spatial  fre¬ 
quencies,  such  t fiat  the  first  harmonic  of  the  grating  transform  is  positioned  at  the  mid¬ 
radius  of  the  detector  ring,  would  yield  a  unique  detection  of  the  grating  frequency  if 
the  sampling  aperture  diameter  and  transform  lens  focal  length  are  properly  selected. 
In  that  section,  we  also  determined  that  the  direction  of  the  Ronchi  ruling  could  he 
determined  uniquely  if  the  grating  was  aligned  perpendicular  to  the  bisector  of  a  de¬ 
tector  wedge.  In  this  section,  we  will  be  more  pragmatic  and  will  consider  the  preliminary 
optical  subsystem  design  of  grating  map  systems.  We  will  see  how  optical  system  para¬ 
meters  affect  the  amount  of  information  that  can  be  recorded,  and  we  will  see  how  the 
interaction  of  parameters  force  the  optical  designers  to  make  trade-offs  for  realistic 
optical  system  specifications.  Yet  we  will  not  have  enough  information  to  design  a  system 
because  of  the  dependence  upon  the  grating  fabrication  capability,  which  will  be  dis¬ 
cussed  in  the  following  section. 


A.  Basic  Optical  Subsystem  Design.  •  Suppose  we  have  a  Ronchi  grating 
map  fabrication  capability  with  a  minimum  writing  spot  diameter  d.  The  maximum 
grating  spatial  frequency  of  the  recorder  system  is  then 

Wm  a  \  =  >/-^1  < 

and  the  minimum  period  of  the  grating  is 

P  =  l/w  =  2d  ’5) 

Knowing  the  dimensions  of  the  frequency  plane  detector,  we  can  use  the  relation 

1  =  R,»ai^,„o 

to  determine  the  transform  lens  focal  length  f  of  the  system,  given  the  mid-radius  of 
the  most  distant  detector  ring  R  and  the  wavelength  X  of  the  laser  light . 


The  minimum  detectable  grating  spatial  frequency  w  jn  is  related  to  the  sampling 
aperture  diameter  I).  To  show  this  relationship,  we  note  that 


w 


m  in 


i  ,v 


m  .i  \ 


(27) 
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where  Pmax  represents  the  Ronchi  grating  period  associated  with  the  minimum  grating 
spatial  frequency  wmjn.  We  know  that  the  Ronchi  grating  period  must  be  less  than  the 
sampling  aperture  diameter  D  to  enable  detection  of  Ronchi  frequency  spectra.  We 
choose 


P„,ax  <  D/N  (38) 

where  the  factor  N  =  2.5  enables  at  least  2.5  periods  of  the  grating  to  be  sampled  with 
safety:  a  factor  smaller  than  2.0  will  result  in  frequency  spectra  of  a  line  ora  line  and 
an  edge.  Thus,  tire  minimum  detectable  Ronchi  grating  spatial  frequency  should  be 


w,»  in  >  2.5/D  (39) 

For  example,  if  1)  =  0.1  mm,  then  a  choice  for  wnijn  would  be  25  lp/mm.  A  clear 
choice  for  wmin  can  be  obtained  by  the  association 


w 


m  in 


(40) 


where  Rmjn  is  the  mid-radius  of  the  detector  ring  closest  to  the  optical  axis  that  will 
satisfy  equation  (39). 


Sampling  aperture  diameter  D  and  transform  lens  focal  length  f  also  influence 
S/N  as  the  simulation  data  in  the  previous  section  indicates.  For  an  f#  less  than  1000. 
i.e. 


f#  =  f/I)  <  1000  (41) 

the  simulation  assures  S/N  >  1.0  lor  WRD-6420  detectors  for  2b  rings  and  32  wedges 

(table  10)  or  27  rings  and  32  wedges  for  the  WRD-6400  detectors  (table  13).  Higher 

... 

1^  combinations  can  be  used  as  indicated  in  tables  13  and  18. 


The  diameter  of  the  sampling  aperture,  coupled  with  detector  responsivity,  will 
determine  the  power  requirements  for  the  system's  laser  light  source.  (Data  are  not 
available  on  laser  power  requirements  at  this  time.)  Available  laser  power,  sample  aper¬ 
ture  diameter,  and  detector  responsivity  will  determine  the  maximum  sample  rates  and 
decision-making  software  that  can  be  optimized  to  operate  within  these  constraints. 
Sampling  and  decision-making  should  be  possible  at  microsecond  rates. 


Selection  of  the  sampling  aperture  diameter  is  of  major  importance  for  the  S/N 
of  the  detected  signal,  and  it  determines  the  resolution  of  the  factor  map  sampling. 
Selection  of  the  sampling  aperture  diameter  should  be  consistent  with  the  precision  of 


the  factor  map.  Table  19  indicates  the  relation  between  selected  sampling  aperture 
diameters  for  various  map  reduction  ratios  and  their  equivalent  diameters  when  projected 
onto  a  map  or  the  ground  surface.  For  example,  a  22"  x  29  (558.8  mm  x  736.6  mm) 
factor  map  having  a  scale  of  1:50,000  when  reduced  lOx  will  have  a  format  size  of 
55.9  by  73.7  mm,  which  is  suitable  for  commercial  70-mm  chips  or  roll  films.  If  a  0.1 
mm  diameter  sampling  aperture  is  used  for  grating  detection  at  this  reduction,  it  has  an 
equivalent  of  50  meters  on  the  ground.  Complete  readout  of  the  factor  map  could  be 
accomplished  with  41  1,612  samples.  Knowing  the  aperture  diameter  and  sample  spacing 
and  knowing  the  format  size  of  the  grating  map  can  be  used  to  establish  specifications 
for  the  optical  or  mechanical  scanning  components. 


B.  System  Design  Trade-Offs.  •  For  any  given  set  of  design  parameters, 
there  is  sufficient  motivation  to  investigate  approaches  that  will  yield  a  more  cost-effec¬ 
tive  system  design.  One  potential  means  for  accomplishing  this  design  is  to  put  more 
factor  maps  on  the  same  recording  material.  For  example,  above  we  indicated  that  a 
70-mm  chip  could  be  used  to  encode  grating  maps  with  lOx  reduction,  and  if  we  now 
consider  35-mm  grating  maps  with  a  20x  reduction,  then  we  could  put  four  grating 
maps  on  a  70-mm  film  chip.  Likewise  16  seperate  16-mm  grating  maps  might  fit  on  a 
70-mm  film  chip  with  40x  reduction,  and  with  65x  reduction  98  seperate  factor  maps 
could  be  encoded  on  a  4  by  5.5-inch  (10  by  14cm)  microfiche  format.  Although 
factor  map  reductions  of  this  magnitude  are  not  practical,  we  will  now  discuss  system 
trade-offs  at  the  designer’s  disposal  that  may  lead  to  smaller  grating  map  formats  for 
selected  applications. 


I .  Focal  Length  Trade-Offs.  Table  19  is  of  value  in  selecting  the  sampling 
aperture  diameter  when  given  a  desirable  level  of  factor  map  reduction  and  the  map 
sample  resolution.  From  the  simulation  data,  we  found  that  the  of  the  optical  system 
(41)  was  a  major  factor  in  determining  the  S/N  for  Ronchi  Grating  Detection.  The  reason, 
of  course,  relates  to  the  aperture  diameter  scaling  of  the  Bessel  function  as  shown  in  table 
1.  Thus,  if  we  achieved  acceptable  S/N  results  for  D  =  0.1  mm  and  f  =  100.0  mm, 
i.e.  f/D  =  1000,  then  we  would  expect  equivalent  S/N  results  for  D  =  0.05  mm 
and  f  =  50.0  mm.  Although  this  is  true,  we  must  also  keep  in  mind  that  given  a  grating 
of  spatial  frequency  w(>  a  change  in  f  will  effect  a  change  in  the  position  of  the  grating 
harmonics  on  the  detector,  e.g. 


r,/f2  =  r,/r2 


(42) 


where  R,  is  the  radius  from  the  optical  axis  to  the  first  harmonic  of  the  grating  function. 
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Or  if  we  fix  the  radius  R0,  such  as  a  mid-ring  radius  of  the  detector,  then  a  change  in 
f  must  be  reflected  in  a  change  in  the  spatial  frequency,  i.e. 

f,  If 2  =  w2/w, 

Thus,  if  we  reduce  the  focal  length  by  a  factor  of  2.0,  we  must  increase  the  grating 
spatial  frequency  by  a  factor  of  2.0  to  keep  the  first  harmonic  of  the  grating  at  radius 
R  We  will  see  in  the  next  chapter  that  this  trade-ofl  is  not  generally  applicable  because 
it  implies  a  capability  for  fabricating  gratings  at  arbitrary  spatial  frequencies. 


2.  Laser  Wavelength  Trade-Offs.  There  is  another  trade-off  available; 
we  can  change  the  wavelength  X  of  the  laser  light  source  within  certain  limits.  As  seen 
from  (36)  and  (40).  the  wavelength  of  the  laser  light  source  is  the  scalar-relating  focal 
length  f,  the  grating  spatial  frequency  w,  and  the  radial  distance  froni  the  optical 
axis  R.  Also,  table  1  indicates  that  the  size  of  the  Bessel  function  for  any  f*  is  directly 
related  to  the  wavelength  X. 


Table  20  lists  some  common  gas  lasers,  their  associated  wavelengths,  and 
a  ratio  value  of  the  wavelength  compared  to  the  common  HeNe  wavelength.  Using  the 
expression  for  ratios  of  the  Bessel  function  in  table  1  and  the  laser  wavelengths  in  table 
20.  we  can  compute  equivalent  sampling  aperture  diameters,  as  shown  in  table  21.  Here 
the  radius  to  the  fourth  Bessel  zeros  are  given  for  detected  values  of  f  and  D  with 
X  =  0.6328  xlO'3  mm.  The  equivalent  D,  for  other  X,,  which  yields  the  same 
Bessel  radii,  are  given  by 

D2  =  D,  (X,/X,)  (44) 

and  are  shown  to  the  right  of  table  21 .  Figure  44  graphs  f/D  versus  radii  from  the  optical 
axis  for  two  selected  Bessel  zeros  and  two  laser  wavelengths.  Changing  from  HeNe  to 
HeCd  lasers  reduces  the  radius  of  the  Bessel  zeros  by  nearly  a  factor  of  2. 


Of  equal  importance,  we  can  improve  S/N  by  selecting  a  shorter  wavelength. 
Modification  of  equation  (41 )  to  reflect  this  sensitivity  to  X  is  given  by 

(X,/X,)  f#  =  (X2/X,)  (f/D)  <  1000  (45) 


When  the  wavelength  of  the  light  souree  is  changed,  it  should  be  remembered 
that  the  mid-ring  spatial  frequencies  will  change  accordingly,  as  shown  in  Table  22 


100 


TABLE  20.  Common  Gas  Lasers  and  Their  Associated  Wavelengths,  X. 


Ratio 

Gas  Lasers  (CM)  Ai/xHeNe 

HeNe  ,  X  =  0.6328  x  10"3  mm  1.000 

Ar  ,  X  =  0.5145  x  10’3  mm  0.813 

1 

x  =  0.4880  x  10“3  mm  0.771 

2 

HeCd  ,  X  =  0.4416  x  10-3  mm  0.698 

l 

x  =  0.3250  x  10-3  mm  0.514 

? 

Gas  Lasers  (Pulsed) 

Ne  ,  X  =  0.5401  x  10‘3  mm  0.854 

N  X  =  0.3371  x  10-3  mm  0.533 

2 
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tor  the  WRD-6420  detector.  Also,  the  change  in  sensitivity  of  t he  ROSA  detectors  at 
shorter  wavelengths  must  be  considered,  which  will  require  higher  laser  power  levels 
for  equivalent  output  signals. 


3.  Trade-Offs  For  Position  Of  Input  Transparency.  Although  the  optical 
configuration  shown  in  figure  4  is  most  commonly  used  in  OPS  systems,  there  is  another 
configuration  of  potential  value.  In  this  configuration,  the  grating  transparency  normal 
to  the  optical  axis  is  placed  behind  the  transform  lens,  as  shown  in  figure  45  (a),  for  the 
grating  positions  2  and  3  instead  of  at  the  conventional  position  1.  Goodman*  discusses 
this  configuration  as  involving  diffraction  of  the  converging  cone  of  rays  in  tne  object 
plane  that  is  detected  and  explained  in  the  focal  plane  by  Fresnel  diffraction.  An  inter¬ 
esting  aspect  of  this  arrangement  is  that  the  aperture  function  transform,  as  indicated 
by  figure  45  (b),  remains  the  same,  no  matter  where  the  transparency  is  placed  between 
the  transform  lens  and  the  detector  and  even  though  the  diameter  of  the  physical  sam¬ 
pling  aperture  decreases  by  a  ratio  of  d/f.  This  aspect  is  due  to  the  constant  f^  of  the 
optical  system.  Moving  the  transparency  behind  the  transform  lens  does  however  change 
the  scale  of  the  spectra  in  proportion  to  d/f.  In  figure  45  (c).  the  first  two  harmonics 
of  a  Ronchi  grating  spectrum  are  shown.  Now  consider  this  grating  to  be  placed  at  one 
of  the  indicated  positions  in  figure  45  (a).  The  physical  location  of  the  harmonics  in  the 
focal  plane  space  is  scaled  according  to 


Thus  in  figure  45  (d,  e,  0.  we  see  the  si/.e  of  the  spectrum  on  the  detector  decreasing 
as  the  transparency  is  moved  closer  to  the  detector.  Also,  the  amplitude  in  the  con¬ 
verging  cone  of  incident  light  will  be  more  intense  by  K.  =  K  f/d  as  the  physical 
sampling  aperture  becomes  smaller.  This  is  also  indicated  in  figure  45  (d,  e,  0. 


Trade-offs  associated  with  this  approach  involve  using  smaller  physical 
sampling  apertures,  with  some  flexibility  of  making  adjustments  or  changing  selections 
versus  the  complexity  of  maintaining  multiple  sets  of  decision  criteria  if  the  frequency 
plane  scaling  flexibility  is  exercised.  Graphically  it  can  be  seen  from  figure  45  that  S/N 
should  decrease  when  figure  45  (b>  is  convolved  with  the  spectra  of  figure  45  (e  and  0 


W .  Goodman.  Introduction  To  Touricr  Optics ,  McGraw-Hill  Hook  Company,  San 
1  ran  cisco  .  ('  A  ,  1968. 
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FIGURE  45.  Gratings  Placed  Before  and  After  Transform  Lenses. 
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for  the  transparency  positions  2  and  3,  which  should  not  be  considered  a  generalization 
until  further  investigated.  Another  possible  limitation  of  this  approach  involves  aber¬ 
rations  in  the  detector  plane  introduced  by  placing  the  transparency  in  the  converging 
beam  of  light.  The  major  aberation  will  be  spherical  aberration,  and  it  should  increase 
in  severity  as  d/f  becomes  smaller. 


Additionally,  the  wavelength  trade-offs  discussed  in  the  previous  section 
can  be  coupled  with  those  of  this  section  to  provide  additional  flexibility  if  desired. 


VI.  FABRICATION  OF  GRATING  MAPS 


In  the  previous  chapter,  the  preliminary  optical  subsystem  design  and  trade-off 
considerations  for  grating  map  systems  were  outlined.  We  saw  that  the  system  depended 
upon  the  grating  map  fabrication  capability.  In  this  chapter,  we  look  at  the  technology 
for  producing  grating  maps  and  complete  the  preliminary  system  design  discussion. 
Then,  approaches  for  extending  capabilities  to  produce  gratings  with  increased  storage 
capacity  will  be  outlined. 


A.  Fabrication  Requirements.  •  For  our  purposes,  a  grating  map  is  an  optical 
transparency  with  factor  map  information  encoded  by  Ronchi  gratings  at  selected  spatial 
frequencies  and  directions.  Producing  a  grating  map  requires  a  capability  to  transform 
factor  map  codings  into  codings  that  can  be  used  to  produce  the  required  gratings.  For 
each  class  of  data  on  the  factor  map,  a  unique  Ronchi  grating  should  be  generated  in 
terms  of  spatial  frequency  and  direction.  Figure  46  indicates  conceptually  how  a  grating 
map  could  be  produced.  Here  the  factor  map  at  the  upper  left  is  transformed  to  a  grating 
map  at  the  lower  right  by  a  technique  that  essentially  involves  spatial  masking  with  a 
seperate  grating  for  each  factor  map  data  class.  In  addition  to  the  grating  map,  there 
must  be  an  associated  legend  indicating  the  correspondence  between  the  gratings  and 
the  factor  map  data  classes.  Also,  a  need  exists  for  a  capability  to  register  automatically 
the  transparency  in  the  optical  system,  e.g.  reference  marks  permitting  automatic  re¬ 
gistration  of  the  grating  map  in  the  optical  system. 


Although  the  grating  map  could  be  produced  by  cutting  masks,  it  would  be  im¬ 
practical  because  (I)  it  is  labor  intensive  and  (2)  it  would  require  a  photographic  re¬ 
duction  capability  to  achieve  a  micrographic  image  for  mass  storage. 
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A  more  practical  approach  would  be  to  use  a  digital  computer  to  program  a  modu¬ 
lated  light  scanning  system  that  would  write  the  gratings  at  the  desired  scale.  We  will 
assume  the  digital  system  can  be  programed  to  accept  factor  map  boundary  data  from 
a  factor  map  and  to  generate  a  signal  for  driving  the  light  modulator  and  scanner  in 
synchronism.  Two  modes  of  scanning  are  available:  vector  scan  and  raster  scan.  The 
vector  mode  requires  a  random  access  capability,  but  provides  faster,  simpler  and  better 
performing  grating  recording  capability:  whereas,  the  raster  mode  is  strictly  sequential, 
but  produces  severe  limitations  on  the  signal  processing  and  writing  speeds. 


Table  23  lists  important  parameters  for  grating  map  recorders.  We  will  briefly 
review  these  parameters. 

Printing  Spot  --  The  size  and  shape  of  the  printing  spot  on  the  recording 
medium  will  largely  determine  the  resolution  of  the  system.  In  general,  the  intensity 
profile  of  the  printing  spot  will  be  Gaussian  and  not  be  a  perfect  square  pulse.  When 
the  beam  shape  departs  from  an  optimum,  it  is  often  necessary  to  provide  a  pixel  overlap. 
The  effective  beam  size  and  shape  is  mainly  determined  by  the  apertures  in  the  system, 
which  act  as  spatial  filters. 


Geometric  Precision  -  Geometric  distortions  of  many  types  plague  scan¬ 
ning  printers.  Deviations  from  orthogonality  of  the  X  and  Y  scanning  axes  cause 
distortions  in  the  output  image.  The  point-to-point  accuracy  of  the  scanner  across  the 
format  determines  the  geometric  precision  of  the  scanner  and  should  be  on  the  order  of 
0.1  of  the  pixel  size.  Rescan  error  is  important  when  two  images  are  superimposed, 
and  line  start  jitter  can  differentially  displace  scan  lines  in  raster  scan  systems. 


Radiometric  Precision  -  The  minimum  intensity  of  the  scanning  system 
determines  the  minimum  density,  or  fog  level,  of  the  image.  This  parameter  becomes 
particularily  important  where  the  rate  of  change  of  intensity  to  distance  is  large.  For 
example,  if  a  grating  with  a  spatial  frequency  approaching  that  associated  with  the  spot 
size  is  to  be  written  in  the  raster  mode  where  the  scan  direction  is  transverse  to  the 
grating,  the  rate  of  change  of  the  intensity  level  will  determine  the  grating  modulation 
at  that  frequency.  If  the  signal  changes  from  'full  on'  to  Tull  off  in  one  pixel  distance, 
the  digital-to-analog  converters  and  the  following  electronics  are  "strained.”  The  band¬ 
width  of  the  system  is  defined  by  the  rate  of  change  of  the  signal  rather  than  the  pulse 
width.  The  radiometric  errors  should  not  exceed  one-half  of  the  least  significant  bit 
of  the  image  data. 


Writing  Speed  -  The  throughput  time  or  writing  speed  to  produce  an 
image  from  start  to  stop  has  been  reduced  significantly  during  the  past  several  years, 
often  at  the  expense  of  geometric  and  radiometric  fidelity. 
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TABLE  23.  Parameters  For  Grating  Recorders. 


Printing  Spot 
Size 
Shape 

Pixel  overlap 

Geometric  Precision 

Orthogonality  of  axes 
Point-to-point  accuracy 
Rescan  error 
Line  start  jitter 

Radiometric  Precision 
Minimum  intensity 
Rate  of  change  of  signal  level 
Radiometric  error 

Writing  Speed 

Operation  and  Maintenance 
Operability 
Maintainabil ity 
Rel iabil ity 
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Operation  and  Maintainencc  -  In  addition  to  being  reliable,  systems  must 
be  simple  to  operate  and  maintain. 


B.  The  State  of  Technology.  •  Table  24  contains  a  partial  list  of  several 
technologies  for  producing  Ronchi  grating  transparencies,  including  examples  of  available 
equipment.  The  completeness  of  the  descriptive  statements  in  terms  of  the  parameters 
of  table  23  illustrate  the  diversity  and  nonuniformity  of  the  specifications  with  which 
manufacturers  describe  their  equipment. 


Although  the  equipment  listed  in  table  24  was  developed  to  generate  hardcopy 
images  from  a  computer,  two  listed  systems  appear  to  have  the  specifications  necessary 
to  make  grating  maps.  The  two  systems  are  Electron  Beam  Recorder  (HBR).  produced 
by  Image  Graphics  Inc.,  and  the  Cathode  Ray  Tube  (CRT)  system,  produced  by  Infor¬ 
mation  International  Inc.  Both  systems  have  raster-  and  vector-scan  modes.  The  high 
resolution  HBR  produces  line  widths  from  3  to  ‘>0  pm  (micrometers)  by  dithering  a  3 -pm 
printing  spot  around  any  of  32,768  x  32.768  addressable  points,  and  the  CRT  system 
adjusts  spot  size  from  3.6  to  19.6  pm  for  vector  scanning.  Both  systems  have  recording 
speed,  format  size,  and  geometric  fidelity  advantages  over  the  other  systems. 


Without  available  test  data,  a  question  arises  of  how  precisely  either  system  could 
produce  a  Ronchi  grating  for  a  selected  spatial  frequency  within  their  operational  range, 
l  or  example,  an  ideal  Ronchi  grating  will  have  sharp  discontinuities  in  its  optical  density 
profiles;  whereas,  gratings  produced  by  Gaussian -scanned  printing  spots  will  have 
smoother  transitions,  and  the  result  will  be  a  degraded  Ronchi  grating.  The  severity  of 
degradation  is  dependent  upon  the  beam  profile.  The  effect  of  this  degradation  is  to 
produce  even  and  odd  harmonics  in  the  grating  spectrum,  which  will  decrease  the  S/N 
of  the  detected  gratings.  However,  a  limited  testing  program  for  generating  gratings  at 
a  few  selected  spatial  frequencies  would  resolve  these  questions. 


C.  Further  Discussion  of  System  Design  Trade-Offs.  •  For  this  preliminary 
report,  we  assume  that  the  above  FBR  and  CRT  systems  will  produce  useful  gratings. 
What  then  can  we  say  about  the  optical  system  for  detection  of  these  gratings  and  how 
many  data  classes  can  we  expect  to  encode  and  decode  with  these  systems?  Although 
it  will  be  a  continuation  of  the  discussions  of  design  trade-offs  (see  section  SB),  this 
section  will  be  concerned  only  with  the  FBR  fabrication  system. 
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TABLE  24.  Examples  of  Grating  Recording  Systems. 
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19.6  um,  64,000  x  64,000  addressable  points,  7.5  cm  x  7.5  cm  display, 
linearity  0.05%,  512  intensity  levels,  0  to  3.15  density  range,  computer 
corrections  for  distortions,  writing  speed  400,000  points/sec. 
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The  minimum  spatial  frequency  grating  generated  by  the  EBR  system  with  a  line 
width  of  90  pm  is  5.5  Line-pairs  per  millimeter  (LP/mm)  and  the  maximum  spatial 
frequency  grating  generated  with  a  3^im-line  width  is  166.7  1/mm.  From  this  infor¬ 
mation,  we  want  to  determine  the  best  system  for  the  parameters  of  sampling  aperture 
diameter  D,  the  focal  length  of  the  transform  lens  f,  and  the  maximum  number  of 
detectable  gratings  corresponding  to  seperate  combinations  of  rings  and  wedges  of  the 
ROSA  detector.  To  assist  this  discussion,  table  25  presents  an  example  worksheet  in 
selecting  the  grating  detection  system  parameter  vrf'W  In  addition  to  restricting  this 
design  to  the  EBR  gratings,  we  assume  that  a  W®.1/  detector  will  be  used  and  that 

a  0.10-mm  sampling  aperture  is  required  t  >  ur-.  '  v.ip:ing  criteria.  However,  the  ap¬ 
proach  is  general,  and  other  criteria  could  be  ft*  u.  -. 


In  Trial  1  of  table  25,  the  objective  «s*  u  an  initial  estimate  for  the  parameter 
values  by  computing  f  from  equation  t3b),  wmjn  from  equation  (39),  and  D  from 
equation  (41).  The  result  is  a  nonstandard  vaJue  for  f,  and  D  was  larger  than  specified. 
In  Trial  2,  a  standarized  value  for  f  (6  inch  «>  is  given,  the  other  values  are  computed, 
and  a  larger  D  value  is  obtained.  In  Trial  3,  values  for  f  and  D  are  given  and  the  value 
f#  is  found  to  be  larger  than  desired.  Reducing  f  to  100.0  mm,  while  keeping 
D  =  0.10  mm  in  Trial  4,  results  in  wmav  at  detector  Ring  32  being  much  greater 
than  the  high  spatial  frequency  EBR  capability,  and  w|njn  for  Ring  7  has  a  spatial 
period  that  violates  equation  (39).  Next.table  2  is  used  to  find  the  detector  ring  having 
a  spatial  frequency  equal  to  or  less  than  the  EBR  high  spatial  frequency  capability,  and 
Trial  5  indicates  this  to  be  Ring  25.  Likewise,  the  lowest  ring  satisfying  equation  (39) 
is  found  to  be  Ring  9  from  table  2.  This  solution  yields  a  total  of  1 7  rings  and  32  wedges 
for  encoding  544  unique  gratings.  The  decoding  system  would  have  a  focal  length  of 
100.0  mm,  and  the  sampling  aperture  would  be  0.10  mm,  as  per  the  specifications. 
The  trials  that  follow  explore  possibilities  of  improving  the  capabilities  by  employing 
different  laser  wavelengths.  As  shown  in  table  25,  Trial  9  achieves  a  similar  number  of 
rings  and  wedges,  but  would  require  a  higher  powered  HeCd  laser  because  of  the  lower 
WRD-6420  sensitivity  at  X  =  0.3250  xlO'3  mm.  Thus,  with  all  trade-offs  considered. 
Trial  5  appears  to  offer  the  best  set  of  preliminary  design  values.  A  similar  approach 
could  be  used  to  investigate  the  applicability  of  the  WRD-6400  detector  or  an  equal 
ring  width  detector  if  the  intended  application  could  support  a  detector  development 
effort.  Likewise,  this  approach  could  be  used  to  design  and/or  assist  in  the  selection  of 
other  viable  systems  from  those  listed  in  table  24. 
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D.  Extending  The  Storage  Capacity  of  Grating  Maps.  •  The  discussion  of 
grating  fabrication  and  the  information  storage  capacity  of  a  grating  map  would  not 
be  complete  without  advancing  some  concepts  that  improve  the  storage  capacity  if 
required.  In  the  preceding  section,  we  developed  an  example  design  showing  that  the 
EBR  system  may  have  the  potential  to  encode  Ronchi  gratings  at  17  spatial  frequencies 
and  32  directions  for  a  total  of  544  thematic  data  classes  for  sampling  with  a  0.10  mm 
aperture.  For  centain  classes  of  data,  this  number  of  data  classes  is  more  than  sufficient, 
but  suppose  the  need  existed  for  producing  grating  maps  with  a  larger  set  of  data  classes. 
This  need  might  arise  if  contour  maps  are  to  be  encoded  for  geographic  areas  having 
severe  relief  display  problems.  Another  need  for  a  large  number  of  data  classes  might 
arise  from  a  requirement  to  combine  factor  maps  so  as  to  keep  the  total  number  of 
graphics  to  a  minimum.  Nevertheless,  these  added  capabilities  will  not  be  gained  without 
paying  some  price.  We  will  look  at  a  few  concepts  for  increasing  capabilities  and  the 
price  to  be  paid  for  the  increase. 


I .  Amplitude  Coding.  All  discussions  have  been  concerned  with  binary 
Ronchi  gratings  in  the  form  of  optical  transparencies,  i.e.  square  wave  gratings  with 
maximum  and  minimum  transmission  values  only.  Table  24  indicates  that  the  EBR 
can  produce  a  modulated  recording  having  64  shades  of  grey  over  an  optical  density 
range  of  0  to  2.3.  Theoretically,  this  method  should  provide  an  encoding  capability  for 
I  7  spatial  frequencies,  32  directions,  and  64  amplitudes  equal  to  1 7  x-32  x  64  =  34,8]  6 
separate  data  classes.  In  reality,  little  margin  is  left  for  noise  associated  with  system 
variability.  Also,  at  low  grating  amplitudes,  very  low  level  OPS  signal  detection  problems 
will  be  experienced.  Thus,  it  is  reasonable  to  select  a  subset  of  the  64  grey  levels,  say  10 
levels,  which  produce  high  S/N  results  similar  to  that  associated  with  the  binary  gratings. 
Just  10  amplitude  levels  may  provide  an  encoding  capability  for  3,481  data  classes  with 
the  EBR.  With  only  two  levels,  the  storage  capability  of  the  binary  gratings  will  be 
doubled. 


To  obtain  additional  storage  capacity  in  this  manner  will  require  a  grating 
preparation  capability  that  assures  the  desired  uniformity.  Also,  the  quality  of  the  signal 
processing  system  must  be  such  that  the  various  signal  levels  can  be  unambiguously 
discriminated.  For  example,  the  laser  light  source  must  have  a  high  degree  of  stability 
over  short  and  long  time  periods.  Thus,  a  simple  laser  for  binary  grating  fabrication 
capability  must  be  replaced  by  a  stabilized  laser  when  the  multi-level  gratings  are  em¬ 
ployed.  Similarly,  the  electronics  must  have  temporal  stability  and  perhaps  more  pre¬ 
cision.  The  software  for  decoding  the  multi-level  grating  detection  data  increases  in 
complexity  and  time  required  per  decision  as  the  number  of  grating  amplitude  levels 
increases. 
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2.  Multiple  Gratings.  Combinations  of  superimposed  gratings,  each  of 
which  is  aligned  so  that  their  separate  spectrums  will  not  fall  on  the  same  detector  wedge, 
provide  another  means  for  increasing  storage  capacity  .  however,  this  approach  is  not  as 
simple  or  powerful  as  the  amplitude  codings.  An  analytic  explanation  of  this  approach 
will  not  be  advanced  here,  but  one  can  easily  conceptualize  that  in  addition  to  the  sep¬ 
arate  spectra  of  the  individual  gratings,  there  will  be  interaction  spectra  that  tend  to  cloud 
the  S/N  mental  picture.  This  approach  could  cause  additional  complexity  in  decision¬ 
making  software  to  determine  whether  superimposed  gratings  or  separate  contiguous 
gratings  existed  in  the  sampling  aperture. 


3.  Color  Gratings.  Amplitude  coding  and  multiple  grating  codes  can  be 
used  with  gratings  prepared  on  color  transparencies.  Separate  lasers  or  narrow  band 
filtered  light  will  be  required  to  decode  the  color  transparency  data.  Table  20  indicates 
a  range  of  blue  through  red  optical  wavelengths  available  in  common  gas  lasers.  Although 
it  would  appear  to  be  promising.this  approach  must  be  tested. 


VII.  SYSTEMS  CONCEPTS 


The  previous  two  chapters  indicated  optical  subsystem  design  and  trade-off  con¬ 
siderations  for  grating  inap  systems,  outlined  the  technology  associated  with  grating  map 
fabrication.  In  this  chapter,  we  outline  more  general  optical  subsystem  concepts  and 
integrate  these  with  the  other  subsystems,  i.e.  the  electronic,  the  digital,  and  the  storage 
and  retrieval  subsystems,  to  present  a  more  unified  view  of  the  major  system  features. 
An  example  system  is  then  described  in  some  detail  to  illustrate  subsystem  interrelations. 


A.  Basic  Subsystems.  •  Important  components  of  the  optical,  electronic, 
digital  and  storage  and  retrieval  subsystems  will  be  reviewed. 

1 .  Optical  subsystem.  The  major  components  of  a  standard  OPS  optical 
subsystem  arc  as  indicated  in  figure  4,  and  the  preceding  two  chapters  addressed  their 
design  considerations  for  OPS  detection  of  grating  maps.  The  basic  optical  configuration 
is  duplicated  in  figure  47(a),  where,  along  with  figure  47(b  and  c).  the  different  con¬ 
figurations  for  placing  a  transparency  in  a  stationary  light  beam  are  indicated.  The  ar¬ 
rangement  shown  in  figure  47(b)  (which  was  given  before  in  V,  B.  3)  and  the  single 
lens  arrangement  shown  in  figure  47(c)  represents  a  further  development  that  is  simpler 
in  concept,  but  will,  in  general,  produce  greater  aberrations  because  of  the  severe  require¬ 
ments  satisfied  only  by  expensive  lens  systems.  One-  or  two-dimensional  sampling  is 
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f(x.y) 


FIGURE  47.  Optical  Arrangements  for  Stationary  Light  Beam 

Ol’S  Systems. 
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accomplished  with  these  optical  configurations  by  translating  the  sampling  aperture 
a(r),  or  more  preferably  by  translating  the  transparency,  while  the  sampling  aperture 
remains  fixed  and  centered  on  the  optical  axis. 


In  general,  stationary  light  beam  OPS  systems,  which  translate  the  trans¬ 
parency  to  obtain  sample  data  over  an  area,  provide  excellent  research  systems,  but 
as  with  the  automated  microdensitometers,  these  arrangements  are  usually  too  slow  for 
operational  applications.  Faster  OPS  sampling  systems  can  be  realized  by  moving  a  tele- 
centric  sampling  beam,  i.e.  a  collimated  beam  parallel  to  the  optical  axis,  over  a  stationary 
transparency.  Two  approaches  will  be  outlined.  The  first  approach,  shown  in  figure 
48(a  and  b)  for  one-  and  two-dimensional  arrangements  respectively,  uses  glass  lenses 
and  galvanometers  to  move  and  form  a  telecentric  scanner,  with  samples  in  plane  1  and 
2.  Here  we  will  briefly  review  the  two-dimensional  arrangement  of  figure  48(b).  The 
beam-forming  optics  (BFO)  are  primarily  responsible  for  producing  a  clean,  collimated 
beam  of  a  given  diameter.  This  beam  is  focused  at  the  center  of  rotation  of  a  galvano¬ 
meter  mirror  MS)  and  this  spot  is  relayed  by  the  relay  optics  (RO)  to  the  center  of 
rotation  of  a  second  galvanometer  mirror  MSl,  which  scans  orthogonal  to  that  of  the 
first.  Generally,  the  galvanometers  are  moved  under  computer  control  so  that  the  sam¬ 
pling  beam  raster  discretely  scans  the  transparency,  and  a  matrix  of  sample  points  are 
obtained;  however,  random  access  scanning  can  also  be  applied.  These  systems  have 
their  limitations  in  that  galvanometer  retrace  time  is  significant,  and  the  galvanometer 
resolution  and  accuracy  must  be  exceptional  to  obtain  an  orthogonal  N  by  N  sample 
matrix  with  a  high  number  of  points  and  with  a  high  reproducibility.  Additionally  if 
the  input  transparency  is  large,  then  the  lenses  must  also  be  large  and  expensive. 


The  second  approach  foi  OPS  sampling  with  sea  inline  telecentric  beams  is 
illustrated  in  figure  49.  Here,  reflecting  optics  in  the  form  of  parabolic  mirrors  reflect 
the  beams  similar  to  the  manner  that  lenses  refract  light.  The  arrangements  in  figure  49 
are  analogous  to  the  three  configurations  using  lenses  in  figure  47.  and  also  to  those 
in  figure  48.  In  figure  49,  the  first  parabolic  mirror  M,  focuses  the  collimated  and 
sized  beam  onto  a  rotation  multi-faceted  mirror  that  reflects  the  beam  to  the  first  mirror, 
where  it  is  collimated,  diffracted  by  the  transparency,  as  in  (a).  Parabolic  mirror  VF 
then  acts  as  a  transforming  element  to  form  the  diffraction  pattern  of  the  sampled 
transparency  on  the  detector  surface.  The  rotating  mirror  is  usually  driven  by  a  stepping 
motor  having  a  high  number  of  resolvable  steps  per  rotation  T.  which  enables  this 
approach  to  have  the  advantage  in  that  it  has  no  retrace  time.  Tilted  facets  on  the  rotating 
wheel  provide  the  scanner  with  a  two-dimensional  raster  capability,  limited  of  course 
by  the  number  of  facets.  Thus,  this  approach  usually  provides  a  better  one-dimensional 
scanner  with  a  motorized  film  drive  mechanism  giving  the  second  and  slower  scan  motion. 
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FIGURE  48.  One-  ami  Two-Dimensional  Telocentric  Beam-Seanning 

Arrangements  Using  Laser. 
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2.  Electronic  Subsystem.  For  our  purposes,  we  separate  the  electrical 
signal  processing  into  electronic  and  digital  subsystems  arbitrarily,  based  upon  whether 
the  signal  is  analog  or  digital,  respectively.  Components  of  the  electronic  subsystem 
relate  principally  to  the  interface  between  the  optical  and  digital  subsystems.  Analytical 
characterization  of  the  ROSA  optical  and  electronic  processing  is  shown  in  figure  5. 
Here  the  detector  changes  the  spatial  light  distribution  in  the  spatial  frequency  plane 
to  an  ordered  set  of  ele  trical  signals,  one  signal  for  each  detector  ring  and  wedge.  These 
signals  must  then  be  amplified  in  a  manner  that  preserves  the  nature  of  the  signal.  The 
signals  can  be  sequentially  amplified,  or  each  signal  can  be  amplified  in  parallel  and  then 
multiplexed  as  shown  in  figure  5.  This  is  followed  by  an  analog  to  digital  (A/D)  con¬ 
verter  that  quantizes  the  analog  detector  element  signals  into  an  ordered  string  of  binary 
signals  that  are  formatted  for  the  digital  computer  system. 


Interfacing  an  OPS  detector  with  a  digital  computer  system  is  not  trivial 
because  the  signal  has  a  potential  dynamic  range  of  106  between  the  high  level  signal 
near  d.c.  and  that  in  the  high  spatial  frequency  areas,  which  is  near  or  within  the  noise 
level  of  the  detector.  The  ROSA  electronic  system  uses  an  autoranging  second  stage 
amplification  that  divides  the  signal  into  decades  as  described  in  Section  3.1  and  more 
specifically  in  equation  ( 10). 


3.  Digital  Subsystem.  The  major  component  of  the  digital  system  is  the 
digital  computer.  The  computer  serves  as  the  heart  and  mind  of  the  system  in  the  sense 
that  it  directs  and  controls  all  operations  within  the  system  and  makes  classification  deci¬ 
sions.  which  are  reflected  in  the  data  or  information  that  are  output  from  the  system.  The 
subsystem  does  only  what  it  is  programed  or  hardwired  to  do;  however,  today’s  systems 
are  sophisticated  relative  to  digital  machines  of  the  not  too  distant  past  and  they  have 
significant  capabilities.  The  digital  subsystems  can  operate  with  today’s  minicomputers. 
Designs  have  flexibility  and  may  choose  to  delegate  selected  processing,  control,  and 
decision-making  responsibilities  to  lower-order  microprocessors  for  one  or  more  reasons. 


Other  components  of  the  digital  subsystem  are  the  input/output  devices 
used  to  communicate  and  interact  with  human  operators. 

4.  Grating  Map  Storage  and  Retrieval  Subsystem.  A  subsystem  concerned 
with  the  storage  and  retrieval  of  the  minified  grating  maps  must  provide  not  only  a 
physical  storage  and  retrieval  capability  but  also  an  index  of  locations  where  grating 
maps  are  stored  and  must  be  returned  to  after  use.  This  indexing  function  can  be  served 
by  data  base  subsystems  that  are  integral  to  the  digital  subsystem. 


Automated  and  semi-automated  components  for  physical  storage  and  re¬ 
trieval  of  large  numbers  of  microfiche  and  smaller  format  transparencies  arc  commercially 
available. 
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B.  Example  System  Concepts.  •  The  above  discussion  of  basic  subsystems 
will  be  supplemented  with  additional  concepts  for  optical  systems.  This  will  be  followed 
by  a  discussion  of  a  simple  system  using  minified  grating  maps,  which  will  integrate  sub¬ 
system  concepts  and  show  how  the  system  is  controlled  and  coordinated. 


1 .  Additional  Optical  Subsystem  Concepts.  Previous  discussions  have 
emphasized  the  potential  for  an  OPS  detector  in  grating  detection.  Within  automated 
systems  context,  it  would  often  be  convenient  to  be  able  to  “see"  the  apertured  portion 
of  the  input  transparency  sampled  by  the  OPS  detector,  e.g.  to  aid  in  resolving  boundaries 
between  contiguous  gratings  within  an  OPS  sampling  aperture,  to  assist  in  registering  the 
image  and  stage  coordinates  prior  to  sampling,  or  to  help  identify  point  and  line  data. 
An  n  x  n  photodiode  array  detector  can  be  added  to  the  OPS  optical  subsystems  in 
figure  47(a)  and  figure  48(a  and  b)  as  shown  in  figure  50.  The  stationary  beam  optical 
subsystem  in  figure  50(a)  is  rather  simple  to  implement,  and  the  focal  length  of  the  in¬ 
verse  transform  lens  is  chosen  so  that  the  sampled  image  fills  the  n  x  n  photodiode 
array  detector.  In  figure  50(b),  a  one-dimensional  telecentric  scanner  may  have  a  gal¬ 
vanometer  scanner  to  direct  properly  the  reconstructed  image  to  the  image  detector, 
and  in  figure  50(c)  the  two-dimensional  telecentric  scanner  requires  two  galvanometers. 
The  corresponding  galvanometer  from  figure  50(b  and  c)  can  be  slaved  to  those  of  figure 
48(a  and  b),  respectively,  to  effect  a  simple  coordinated  scan. 


Figure  51  indicates  how  multiple  factor  maps  of  the  same  geographic  area 
may  be  sampled  in  a  coordinated  fashion.  This  is  a  simple  extension  of  the  system  shown 
in  figure  47(a). 


Two  example  illustrations  of  how  a  spatial  light  modulator  (SLM)  may  be 
used  in  an  OPS  system  are  shown  in  figures  52  and  53.  Several  types  of  SLM  are  currently 
under  study  at  a  number  of  laboratories  for  applications  in  coherent  optical  image  pro¬ 
cessing  systems.  Casasent6  presents  a  review  of  the  technology.  For  the  applications  in 
these  illustrations,  an  SLM  such  as  the  ltek  Corporation  Pockel's  Readout  Optical 
Modulator  (PROM)  device  with  record,  read,  and  erase  capabilities  and  properties  of  high 
dynamic  range  and  analog  image  manipulation  with  bias  adjustments  can  be  used. 


6D.  Casasent,  “Spatial  Light  Modulators,"  Proc.  IL'KE,  Vol.  65,  No.  1,  Jan  1  977  , 
pp.  14  3-1  5  7. 
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FIGURE  50.  Addition  of  a  Photodiode  Array  Detector  to  Enable 
“Seeing”  of  the  Image  Sampled  by  the  OPS  Detector. 
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FROM  BS1  Sl[l  L, 

laser  ’ 


FROM  fiS4  S4 

LASER  L, 


SIDE  VIEW  SCHEMATIC  IN  X  -  DIRECTION 


DW/TS 


i  =  1.  2,  3.  4 

BS,  -  Beam  Shaping  Optics 
S,  -  Shutters 
TS  -  Translation  Stage 
Lj  -  Transform  Lens 
Mj  -  Mirror 
SM  -  Stepping  Motor 
Mj  -  Spinning  Mirror 
Lj  -  Imaging  Lens 
D  -  OPS  DETECTOR 


operational  sequence 

Grating  to  be  sampled  is  selected  by  Computer  (example,  I  |,  then 
Shutter  S(  is  opened  and  S2.  S3,  and  S_,  are  closed.  Simultaneous!.  , 
the  stepping  motor  rotates  Mirror  Ms  into  position  so  that  lens  L|  car. 
image  the  OPS  from  I  onto  the  Detector,  D.  Detector  D  is  mounted  on 
a  translation  stage  that  is  serve  driven  to  center  tag  Detector  on  the  optical 
axis  and  thus  compensate  for  alignment  errors  in  positioning  the  rotating 
Mirror,  M,. 

The  registered  gratings  are  mounted  on  a  translation  stage  which 
moves  orthogonally  so  as  to  allow  four  sample  beams  to  illuminate  the  same 
grating  coordinates. 


FIGURE  51,  Concepts  for  Multiple-Coordinated  Factor  Map  Gratings. 
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FIGURE  52.  Example  of  an  OPS  Configuration  Using  a  Spatial  Light 
Modulator,  a  Rotating  K  Mirror,  and  a  One-Dimensional 
Photodiode  Array  Detector. 
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FIGURE  53.  I  xample  of  an  OPS  Configuration  Using  a  Spatial  Light 
Modulator,  a  Two-Dimensional  Scanner  and  a  Two-Dimen¬ 
sional  Photodiode  Arrav  Detector. 


The  example  shown  in  figure  52  has  an  SLM  that  records  the  spectrum 
of  the  sampled  portion  of  the  transparency.  The  SLM  is  a  reflection  device  with  a 
dichroic  filter  on  its  face  and  is  nearest  to  the  transform  lens.  This  filter  passes  the  OPS 
record  light  of  wavelength  Aj ,  while  serving  as  a  reflector  for  the  collimated  readout 
light  A.,.  The  light  is  doubly  modulated  image-wise  by  the  OPS  record  during  its  two 
passes  through  the  SLM  crystal.  A  rotating  K  mirror  serves  to  rotate  the  spectrum  on 
a  one-dimensional  detector.  Decision  criteria  associated  with  this  system  would  determine 
the  spatial  radius  and  direction  of  the  first  harmonic  of  the  grating  being  sampled.  The 
K  mirror  rotates  so  as  to  keep  the  optical  axis  stationary,  with  the  image  rotating  at 
twice  the  angular  velocity  as  the  prism. 


For  the  example  shown  in  figure  53,  the  OPS  pattern  is  recorded  in  the  SLM. 
By  changing  the  electrical  bias  on  the  device,  the  positive  pattern  is  changed  to  a  negative 
pattern  and  thresholded.  A  scanning  spot  traces  a  high-resolution,  two-dimensional  raster 
over  the  SLM  and  light  is  transmitted  through  the  device  only  above  the  negative  thresh¬ 
old  so  that  detection  and  identification  of  the  grating  direction  and  spatial  frequency 
is  rather  simple  with  a  two-dimensional  photodiode  array  detector.  In  this  illustration, 
a  transmission  SLM  is  shown,  but  a  reflection  SLM  could  be  used  by  moving  all  com¬ 
ponents  associated  with  the  scanner  to  the  opposite  side  of  the  SLM  as  in  figure  52. 


2.  A  Simple  System  Example.  A  simple  automated  system  for  machine 
readable  factor  maps  is  shown  in  figure  54.  The  purpose  of  this  discussion  is  to  indicate 
the  coordination  of  the  previously  discussed  components  of  an  automated  OPSA  system 
and  to  present  some  additional  concepts  not  developed  earlier.  We  will  first  briefly  review 
the  system  components  and  how  they  are  coordinated,  then  we  will  present  a  scenario 
related  to  the  possible  use  of  such  a  system. 


Most  components  of  the  optical  subsystem  in  figure  54  have  been  discussed 
previously  except  for  the  variable  optical  density  wedge  (2)  and  laser  intensity  level 
detector  (4),  which  may  serve  in  combination  to  monitor  and  control  the  laser  output 
by  feedback  through  the  computer.  The  rotatable  disc  with  sampling  apertures  (6) 
also  operates  under  computer  control,  as  does  the  two-dimensional  translation  stage  (7). 
The  grating  map  storage  and  retrieval  subsystem  (8)  is  similar  to  a  carousel  that  contains 
grating  map  transparencies  with  more  than  one  grating  map  per  transparency.  This  sub¬ 
system  also  has  capability  for  loading  and  unloading  the  translation  stage  with  trans¬ 
parencies  under  computer  control. 


Electronic  subsystem  components,  except  the  detectors,  have  not  been 
shown  in  figure  54  so  as  to  keep  the  illustration  simple  and  uncluttered.  The  two  de- 
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SYSTEM  COMPONENTS 

1.  LASER 

2.  Variable  Optical 

Density  Wedge 

3.  Beam  Splitter 

4.  LASER  Intensity 

Level  Detector 

5.  Beam  Forming  Optics 

6.  Rotatable  Disc  with 

Sampling  Apertures 

7.  Grating  Map  Transparency 

on  Two-Dimensional 
Translation  Stage 

8.  Grating  Map  Storeage 

and  Retrieval  System 

9.  Transform  Lens 

10.  Beam  Splitter 

11.  OPS  Detector 

12.  Inverse  Transform  Lens 

13.  Photodiode  Array 

Detector 

14.  Micro-Computer  for  1 1 

15.  Micro-Computer  for  13 

16.  Mini-Computer 

17.  Memory  Devices 

18.  Output  Devices 

19.  Input  Terminal 


FIGURE  54.  A  Simple  Automated  OPSA  System  for  Reading 
Optically  Stored  MGI  Factor  Maps. 
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tectors.  similar  to  those  in  figure  50(a).  are  used  for  sampling  the  optical  power  spectrum 
with  a  ROSA  detector  (11)  and  a  two-dimensional  photodiode  array  detector  (13)  for 
sampling  the  input  image. 


A  microcomputer  (14  and  15)  can  be  interfaced  directly  to  each  detector 
to  improve  system  efficiency  by  removing  processing,  control,  and  decision  tasks  from 
the  central  computer.  These  microcomputers  can  be  micro  programed  to  accomplish 
the  data  acquisition  and  formatting  temporary  storage  in  local  buffers,  the  low  level 
data  processing  and  decision-making,  and  the  data  transfer  directly  to  memory  via  a  DMA 
(direct  memory  access).  The  microcomputers  are  controlled  by  a  central  minicomputer 
(16)  that  serves  to  coordinate  and  control  subsystem  operations  and  to  process  data  for 
internal  use  or  to  satisfy  requirements  placed  in  the  system.  The  central  computer  con¬ 
trols  the  laser  output  level  by  periodically  sampling  the  laser  intensity  level  detector  (4) 
signal  after  A/D  electronic  conversion  (A/D  converter  not  shown).  The  central  computer 
then  translates  the  variable  optical  density  wedge  via  a  solenoid  or  stepping  motor  to 
effect  the  desired  light  level  input  for  the  system.  It  also  controls  the  movement  of  the 
two-dimensional  translation  stage  (7).  the  rotatable  disc  with  sampling  apertures  (6). 
and  the  grating  map  transparency  storage  and  retrieval  system  (8).  The  central  computer 
accepts  a  problem  statement  for  a  geographic  area  (19)  and  selects  appropriate  factor 
map  transparencies  with  the  aid  of  a  data  base  management  system  and  other  collateral 
data  in  memory.  The  selected  solution  approach  will  dictate  the  geographic  areas  to  be 
scanned,  the  sample  spacing,  and  the  sampling  aperture  diameter.  The  central  computer 
then  controls  movement  of  the  two-dimensional  translation  stages  (7).  the  rotatable 
disc  with  sampling  apertures  (6).  and  the  grating  map  transparency  storage  and  retrieval 
system. 


The  spatial  image  detector  is  used  to  register  the  transparency  in  the  stage 
coordinate  system.  The  grating  legend,  stored  in  memory  for  the  sampled  transparency, 
provides  information  to  decide  grating  spatial  frequency  and  direction  to  obtain  data 
classifications.  These  decisions  are  processed  in  MGI  application  programs  and  the  results 
are  output  in  readable  form  ( 1 8). 


A  short  scenario  of  the  manner  in  which  the  system  of  figure  54  may  be  used 
in  a  MGI  application  problem  follows.  The  problem  statement  is  presented  to  the  system 
by  an  operator  through  the  input  station  (19)  by  problem  type  with  selected  parameters 
defined,  e.g.  for  a  cross-country  mobility  problem,  the  parameters  may  be  initial  geo¬ 
graphic  location,  final  destination,  type  of  vehicle,  number  of  vehicles,  etc.  Application 
software  then  selects  types  of  data  desired  and  checks  the  data  base  to  determine  if 
these  data  types  are  available  for  the  problem  area. 


129 


Suppose  this  concerned  MGI  Terrain  Factor  Maps  stored  as  Ronchi  grating 
maps  on  transparencies  is  an  automated  storage  and  retrieval  system  as  in  figure  54. 
The  computer  system  then  controls  retrieval  of  a  selected  factor  map  (transparency) 
from  the  storage  carousel  and  inserts  it  into  a  holder  on  the  two-dimensional  translation 
stages.  A  short  automated  registration  procedure  is  accomplished  with  the  aid  of  the 
photo  diode  array  detector  to  align  the  transparency  coordinate  system  with  the  trans¬ 
lation  stage  coordinate  system.  Now  the  stage  coordinates  can  be  computed  to  obtain 
the  desired  transparency  sampling.  Each  time  the  stages  are  incremented  to  a  new  sam¬ 
pling  location,  the  ROSA  detector’s  microcomputer  is  directed  to  obtain  a  set  of  de¬ 
tector  element  samples.  While  the  stage  is  translating  to  the  next  sample  location,  the 
sampled  data  can  be  decoded  and  only  the  factor  map  class  symbol  is  stored  temporarily 
in  a  buffer  within  the  microcomputer.  These  processed  data  are  then  block  transmitted 
to  central  memory  via  direct  memory  access  when  the  buffer  is  full  or  the  sampling  com¬ 
pleted.  Other  factor  map  data  for  the  same  area,  e.g.  soils,  slopes,  vegetation,  etc.,  can  oe 
processed  in  like  manner.  The  central  computer  can  then  use  these  spatial  data,  along 
with  temporal  data,  e.g.  local  weather  history  and  projected  weather,  to  compute  or 
decide  on  the  cross-country  mobility  data  class  for  each  sample  location.  Processed  data 
can  then  be  formatted  and  output  on  a  monitor  and/or  plotter  for  inspection  and/or 
distribution. 


A  somewhat  similar  program  to  the  system  in  figure  54  is  flow  charted  in 
figure  55.  This  simple  example  relates  to  only  one  of  many  possible  ways  in  which  an 
automated  grating  map  storage  and  retrieval  capability  may  be  employed  as  a  machine- 
readable  MGI  factor  map  mass  store  system.  These  concepts  must  now  be  tested  and 
verified  for  correctness  as  well  as  completeness.  The  potentials  at  this  preliminary  stage 
are  quite  promising. 


VIII.  SUMMARY 


This  report  analysized  optical  power  spectral  (OPS)  approaches  to  machine- 
readable  factor  maps  for  use  in  terrain  information  systems.  The  objective  was  to  offer 
a  more  compact,  simpler,  and  cheaper  factor  map  data  base  storage  and  retrieval  systems 
than  would  be  available  by  all  digital  technology.  The  basic  concept  of  the  system  in¬ 
volves  use  of  OPSA  for  decoding  data  from  reduced  factor  maps  that  have  been  encoded 
by  spatial  frequency  and  direction  of  Ronchi  optical  gratings. 


130 


Revise  Solution 

Stragety 

? 

Yes 

No 

1 

iH 

' 

Transmit 
Buffer  Con¬ 
tents  to 
Memory  via 
DMA 

1 

Store 

In 

Buffer 

FIGURE  55.  Program  Associated  With  Automated  OPSA  System 

in  Figure  54. 
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Aii  OPSA  system  consists  of  a  coherent  optical  sampling  subsystem,  an  electronic 
detection  and  signal  processing  system,  and  a  digital  computer  with  associated  software 
for  control  and  decision-making.  In  Chapter  111,  an  analytical  characterization  of  an 
OPSA  system  was  discussed  when  the  input  transparencies  are  ideal  Ronchi  optical 
gratings.  Particular  attention  was  given  to  characterizing  the  resultant  optical  signals 
from  ROSA  (Recording  Optical  Spectrum  Analyzer)  photodiode  detectors  having  32 
concentric  ring  detectors  on  one  semi-circle  and  32  equally  sized  wedge  detectors  on  the 
conjugate  semi-circle.  This  result  formed  the  basis  for  presentating  a  computer  simulation 
program  in  the  following  chapter,  along  with  simulation  results  from  its  use  under  a  wide 
range  of  operational  parameters  with  two  commercially  available  ROSA  detectors  of 
slightly  different  detector  element  configuration.  One  detector  had  mid-radii  of  the  ring 
elements  spaced  logrithmically.  and  the  other  detector  had  linear  spacings.  These  simu¬ 
lation  results  indicated  that  as  many  as  27  rings  and  32  wedges  yield  factor  map  encoding 
potentials  for  864  thematic  data  classes  (27  x  32)  with  detection  signal-to-noise  ratios 
greater  than  1 .0  when  the  f#.  i.e.  the  ratio  of  the  transform  lens  focal  length  to  the 
sampling  aperture  diameter,  is  less  than  1,000.  These  results  motivated  an  attempt  to 
redesign  the  ROSA  detector  element  configuration  to  effect  better  performance.  The 
detector  redesigns  witli  equal  ring  element  spacings  resulted  in  a  capability  of  detecting 
30  rings  and  32  wedges,  i.e.  660  thematic  data  classes,  with  a  signal-to-noise  ratio  greater 
than  1 .0. 


Preliminary  design  of  an  optical  subsystem  for  detection  of  grating  maps  is  de¬ 
veloped  in  Chapter  V  to  indicate  how  the  interaction  of  system  parameters  forces  the 
designer  to  make  trade-offs  for  realistic  optical  system  specification.  Here  particular 
attention  is  given  to  sampling  aperture  diameter,  which  determines  the  resolution  of 
factor  map  sampling.  For  example,  a  22-  by  26-inch  factor  map  having  a  scale  of 
1:50.000  when  reduced  lOx  will  have  a  format  size  of  55.9  mm  by  73.7  mm,  which 
is  suitable  for  commercial  70-mm  chips  or  roll  films,  and  a  0.1-mm-diameter  sampling 
aperture  will  have  an  equivalent  resolution  of  50  meters  on  the  ground.  Trade-offs  leading 
to  smaller  grating  map  formats  are  then  investigated.  Significant  findings  indicate  that 
a  0.05 1 -mm  sampling  aperture  having  approximately  a  25-meter  ground  resolution  can 
have  the  same  signal-to-noise  results  as  the  0.1  mm  aperture  with  50-meter  ground  re¬ 
solution  if  the  laser  light  wavelength  is  changed  from  the  HeNe  0.6328pm  to  HeOd 
at  0.350ptn. 


In  Chapter  VI,  the  technology  for  producing  grating  maps  is  reviewed  and  ap¬ 
proaches  for  extending  capabilities  to  produce  gratings  with  increased  storage  capacity 
are  outlined.  Two  systems  appear  to  have  the  specifications  necessary  for  generating 
grating  maps:  the  electron  beam  recorder  made  by  Image  Graphics  Inc.,  and  the  cathode 
ray  tube  system  available  from  Information  International  Inc.  However,  without  test 
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data  there  is  a  question  of  how  precisely  either  system  could  produce  a  Ronchi  grating 
for  a  selected  spatial  frequency  within  their  operational  range.  An  example  worksheet 
(table  25)  is  given  to  outline  selection  of  the  grating  detection  system  parameter  values, 
which  indicates  that  the  election  beam  recorder  has  a  potential  to  encode  Ronchi  gratings 
at  17  selected  spatial  frequencies  and  32  directions  for  a  total  of  544  thematic  data 
classes.  Means  for  increasing  the  number  of  thematic  data  classes  are  discussed  for  am¬ 
plitude  coding,  e.g.  with  10  gray  scale  levels  for  the  Ronchi  gratings  the  electron  beam 
recorder  could  encode  3,481  separate  thematic  data  classes  in  a  single  grating  map.  Like¬ 
wise.  multiple  gratings  composed  of  superimposed  gratings  or  gratings  prepared  in  color 
transparencies  decoded  by  selected  laser  light  wavelength  illumination  provide  potential 
for  increasing  data  capacity  in  a  grating  map. 


In  Chapter  VII,  concepts  for  the  optical,  electronic,  digital,  and  storage  and  re¬ 
trieval  subsystems  are  described.  For  example,  several  optical  subsystem  configurations 
are  described  and  illustrated  for  ( 1 )  stationary  light  beam  optical  subsystems  with 
moving  input  transparencies  and  (2)  a  telecentric  scanning  beam  optical  subsystem  with 
stationary  input  transparencies.  It  is  shown  that  one-or  two-dimensional  telecentric 
scanning  can  be  accomplished  with  either  the  standard  refractive  optics  afforded  by 
glass  lenses  or  by  reflective  optics  with  parabolic  mirrors.  Integration  of  a  photodiode 
array  detector  into  an  OPSA  optical  subsystem  is  presented  to  assist  in  registering  image 
and  stage  coordinate  systems  as  well  as  serving  as  a  possible  aid  in  resolving  boundaries 
between  contiguous  gratings  on  an  input  grating  map.  It  is  shown  how  spatial  light 
modulations,  such  as  Itek's  PROM  device,  might  be  used  in  OPSA  systems. 


Also  in  this  Chapter,  coordination  of  the  various  subsystems  into  an  integrated 
automated  system  is  described  in  an  example  scenario.  This  system  features  a  carousel 
for  storing  the  reduced  factor  map  data  encoded  as  Ronchi  gratings.  The  simple  system 
uses  a  minicomputer  as  a  central  processor  to  coordinate  and  control  the  various  sub¬ 
system  operations  and  to  process  data  for  internal  use.  Microcomputers  are  interfaced 
to  an  OPS  detector  and  to  an  image  plane  photodiode  array  detector  to  increase  system 
efficiency  by  removing  processing,  control,  and  decision  tasks  from  the  central  computer. 
OPSA  samples  can  be  preprocessed  and  classified  within  the  microprocessors,  and  the 
reduced  class  ncmonic  can  be  stored  in  local  buffers  until  transmission  to  central  system 
memory  via  direct  memory  access  when  the  buffer  is  full.  This  system  could  be  pro¬ 
gramed  to  operate  in  a  true  automated  fashion  by  accepting  formatted  problem  state¬ 
ments  for  a  defined  geographic  area  and  for  a  defined  problem  for  which  an  application 
algorithm  exists,  e.g.  cross-country  mobility  between  selected  points  for  a  given  class 
of  vehicles.  The  central  computer  could  select  the  appropriate  factor  map  transparencies, 
transfer  them  sequentially  to  the  sampling  system,  register  them  in  the  sampling  coor¬ 
dinate  system,  and  select  the  areas  of  the  transparencies  to  be  sampled,  the  sample 
spacing,  and  sampling  aperture  size.  The  OPSA  detector  microcomputer  could  process 
sample  data  while  the  translation  stage  is  moving  to  the  next  sample  position.  The  central 
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computer  could  use  these  spatial  data,  along  with  temporal  data,  e.g.  local  weather 
history  and  projected  weather,  to  compute  or  decide  on  a  problem  answer  for  each 
sample  location.  The  processed  data  could  then  be  formatted  and  displayed  for  human 
inspection  or  placed  in  a  computer  file  for  access  for  map  overprinting  by  automated 
cartography  systems. 


IX.  CONCLUSIONS 


This  preliminary  investigation  of  optical  power  spectral  analysis  approaches  to 
machine-readable  factor  maps  for  use  in  terrain  information  systems  has  developed 
information  and  data  indicating  that  the  concepts  are  feasible,  experimental  testing  of 
the  simulated  system  concepts  can  provide  data  useful  for  system  design  and  imple¬ 
mentation  planning. 


134 


